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A thorough study and error analysis was made of the Van den Akker or "auxiliary sphere" method of 
determining a scale of directional-hemispherical reflectance factor. The effects of a non-Lambcrtian distribution of 
the reflected radiation, including retroreflection, were included in this study. Three working standards were 
measured to an uncertainty in reflectance of less than ±0.0015 and these will be used as a basis for a new, more 
accurate NBS scale of 6°-nemisphericaI reflectance factor. The new scale and the NBS scale established in 1965 
are in agreement to within the uncertainty of ±0.005 assigned to the 1965 scale. 
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Foreword 



We have arranged this paper in a way whieh should 
accomodate readers who have an interest in this work from 
three quite different points of view: 

1. Those who must make decisions based upon the results 
of the measurements are addressed principally in the 
Summary which directly follows this foreword. The 
material in parts I. and V.B. would also be of interest 
to such readers. 



Applications in which the actual value of the reflec- 
tance is of secondary importance, but for which the 
basis of measurement must be extremely stable in time. 
These applications include quality control in auto- 
mated production and specifications involving color or 
appearance of finished products. For such applica- 
tions, the instruments which are used are calibrated 
with material standards, and the laboratories supplying 
these standards must be able to measure reflectance 
directly in order to measure the standards and verify 
their stability. 



2. Those who wish to use the Van den Akker auxiliary 
sphere method for determining absolute reflectance 
should find the material in parts II, III, and V.A. to be 
of particular interest, and part IV can be read superfi- 
cially, if at all, in order to glean a few technological 
ideas. 

3. Those who are interested in a detailed understanding of 
the way in which these measurements are carried out at 
NBS and in a detailed discussion of the error analysis 
should, after a careful reading of parts II and III, place 
their emphasis on part IV. 

Summary 

Of the light or other optical radiation incident upon a 
surface in a given direction and at a given wavelength, a 
fraction is reflected from the surface. This fraction is called 
the spectral directional-hemispherical reflectance of the sur- 
face, and measuring it accurately is important in two different 
classes of applications: 

1. Applications in which the value of the reflectance is of 
direct importance. Such applications include radiative 
energy transfer as in solar energy devices, lighting 
system engineering, calibrating radiometers in remote 
sensing satellites, and formulating the pigmentation in 
paints and other finishes. 



The economic benefits derived from having this portion of the 
measurement system under control are very large, coming in 
the form of an accumulation of modest benefits over a very 
large base of application [l]. 1 In order to realize these 
advantages, it is important that the uncertainty in the mea- 
surements be commensurate with the uniformity and stability 
of the reflectance of the surfaces encountered in practice. 

For a number of applications, it is important to measure 
directional-hemispherical reflectance factors near 1 with an 
uncertainty in the neighborhood of ±0.001. The present 
extensive work on spectral directional-hemispherical reflec- 
tance was undertaken at NBS for two reasons. First, the 
stated uncertainty of ±0.005 for the NBS reflectance factor 
scale established in 1965 [2] 1 was too large for many of the 
applications for which we were called upon to standardize the 
measurements. Second, and even more disturbing, intercom- 
parisons between the scales of national standardizing labora- 
tories in connection with the work of the International Stand- 
ardization Organization (ISO) revealed differences as great as 
0.015 between the measured value of reflectance of the same 
samples. Such a large difference can have serious economic 
consequences in international trade in finished goods such as 
paper. 

As a first step in this work, we have investigated in great 
detail the Van den Akker auxiliary sphere method of deter- 
mining spectral directional-hemispherical reflectance which 



1 Figures in brackets indicate the literature references on page 49. 
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has been used in the past at NBS, improving upon the 
techniques and providing the necessary corrections to make 
the measurements precise and accurate to within ±0.0015 to 
the best of our knowledge. This step has been completed and 
is reported in this paper. The new and former NBS scales of 
measurement agree to well within the combined measurement 
uncertainties. 

As a second step, we plan to investigate the techniques 
used by the other major national laboratories and to work with 
our colleagues in these laboratories to determine the cause of 
the discrepencies which have been encountered internation- 
ally. As part of this step, we have already completed a set of 
measurements using a second method of determination, the 
Sharp-Little method. The results obtained with that method 
are in good agreement with those reported in this paper and 
are to be published soon in another paper in this same series. 
We have visited the National Research Council Laboratories 
in Canada (NRC) for detailed discussions of these results. We 
also plan to conduct experiments with the Korte method 
currently used by the Physikalisch-Technische Bundesan- 
staldt (PTB) in Germany and possibly with one or two other 
approaches less commonly used. 

The scale of directional-hemispherical reflectance as cur- 
rently established will be disseminated through standards 
supplied through the NBS Office of Standard Reference 
Materials and through commercial secondary standards labo- 
ratories in the United States. The improvements in tech- 
niques which have been developed as a result of this work 
will be submitted to the appropriate committees of ASTM, 
ANSI, and TAPPI for possible incorporation in standard 
procedures. When the international discrepancies have been 
eliminated, we will be working with NRC, PTB, and possibly 
other laboratories as standardizing laboratories for ISO re- 
flectance measurements. 

I. Introduction 

Directional-hemispherical (d/h) reflectance factor mea- 
surements are important in a wide variety of applications. If 
the results of these measurements are to be a useful tool for 
technical communication, the measurements must be made 
accurately. Most reflectometers are not capable of measuring 
d/h reflectance factor directly, but can only compare the 
reflectance factors of two objects. The calibration of such 
instruments is accomplished by measuring a standard object 
which has a known reflectance factor. To see that accurately 
measured reflectance standards are available to the measure- 
ment community is one of the primary responsibilities of the 
spectrophotometry group of the Radiometric Physics Section 
of the Institute for Basic Standards. 

In the development and production quality control of fin- 
ished products in which appearance is an important factor, 
the measurement of reflectance should be accurate to within 
±0.002. This level of accuracy or better is also important to 
rapidly evaluating the stability of reflecting materials under 
weathering and ageing. These two types of applications are 
the ones which commonly call for the lowest measurement 
uncertainty. It is difficult to produce highly reflecting sur- 
faces for which the reflectance is reproduced to better than 
±0.001 and the reflectance of most surfaces is not even 
uniform to this degree. Therefore, a reflectance measuring 
capability for which the uncertainty is less than ±0.001 is 
both necessary and sufficient for a national standardizing 



laboratory. The work described in this technical note is part 
of an effort to reduce the uncertainty in diffuse reflectance 
factor measurements at NBS from an estimated ±0.005, 
which it has been in the recent past, to ±0.001. 

The National Bureau of Standards (NBS) has established 
its scale of d/h reflectance factor in 1965 through an exten- 
sive series of measurements by Goebel, Caldwell, and Ham- 
mond [2]. At that time, Vitrolite [3] glass standards [4] to 
calibrate the General Electric Recording Spectrophotometer 
(GERS) [5] (Cat. 5962004 G28 No. 732986) were measured. 
Until recently, that instrument has been used for most reflec- 
tance measurements made at NBS. The reflectance of the 
Vitrolite standards has been shown to be very stable by 
measurements made over a period of thirty years relative to 
freshly prepared MgO surfaces [6]. In 1974, the scale of 
measurement was rechecked using the same apparatus used 
in the 1965 experiments and the agreement was within the 
experimental error associated with the measurements. 

In April of 1974 and through the following year, it was 
determined from measurements made on a number of samples 
that there was a systematic difference between the scales of 
measurement being used by NBS and the National Research 
Council Laboratories of Canada (NRC). This difference was 
approximately 0.015 at the short wavelength end of the 
visible spectrum and decreased more or less regularly to 
approximately 0.01 at the long wavelength end of the spec- 
trum (table I). A similar intercomparison between NRC and 
Physikalisch-Technische Bundesanstalt of Germany (PTB) 
[7] revealed only slight differences between the measure- 
ments being made by these laboratories (table II). Since the 
NBS uncertainty at that time is conservatively estimated to be 
±0.005 and the NRC uncertainty is conservatively estimated 
to be ±0.003, the difference is clearly significant. 

Table I 

Data from an intercomparison between NRC and NBS of reflectance 
measurements on a sprayed BaSC>4 coating*. (June 1975) (Spectral direc- 
tional (6°)-hemispherical reflectance). 







Reflectance 




Wavelength 


NRC 


NBSf 


Difference 
(NRC-NBS) 


400 


0.961 


0.975 


0.014 


420 


.963 


.975 


.012 


440 


.966 


.978 


.012 


460 


.970 


.980 


.010 


480 


.970 


.981 


.011 


500 


.972 


.981 


.012 


520 


.972 


.985 


.013 


540 


.972 


.985 


.013 


560 


.974 


.984 


.010 


580 


.975 


.985 


.010 


600 


.975 


.985 


.010 


620 


.976 


.982 


.006 


640 


.975 


.986 


.011 


660 


.976 


.985 


.009 


680 


.976 


.985 


.009 


700 


.977 


.985 


.008 


720 


.978 


.985 


.008 


740 


.979 


.985 


.007 



* Samples prepared at NRC. 

t Data from NBS test 232.14/49D. 



The methods used by NBS, NRC, and PTB are all differ- 
ent. In any such case of disagreement between measure- 
ments, all measurements and methods are logically suspect 



30 



Table II 

Data from an intereomparison between NRC and PTB of reflectance 
measurements on a pressed BaS04 tablet [7] (Spectral Hemispherical- 
directional (0°) reflectance factor). 



Reflectance 


Wavelength 


NRC 


PTB 


Difference 
(PTB-NRC) 


370 


0.961 


0.965 


0.004 


380 


.969 


.971 


.002 


390 


.974 


.976 


.002 


400 


.977 


.979 


.002 


420 


.982 


.983 


.001 


440 


.985 


.984 


-0.001 


460 


.986 


.985 


-0.001 


480 


.987 


.986 


-0.001 


500 


.987 


.987 





550 


.988 


.988 





600 


.988 


.988 





650 


.987 


.988 


0.001 


700 


.987 


.989 


.002 


750 


.987 


.988 


.001 



until the cause of the difference is located. However, since 
the NRC and PTB measurements were in close agreement, it 
seemed reasonable to investigate the method used by NBS 
first. A Diffuse Transmittance and Reflectance Reference 
Spectrophotometer [8], which will be referred to by the 
acronym DRS for convenience, was completed at NBS in 
October 1975. This has been used to make detailed measure- 
ments of all aspects of the Van den Akker auxiliary sphere 
method for realizing an absolute scale of d/h reflectance 
factor [9], which was the method used to establish the NBS 
scale. The results of this investigation are reported in this 
paper. 

As a result of this investigation a new NBS scale of diffuse 
reflectance factor measurements for 6° incidence and hemi- 
spherical collection has been established which is believed to 
be accurate to within ±0.0015. The new scale of reflectance 
factor differs by less than 0.002 from the scale formerly used 
by NBS. Since this investigation revealed no cause for the 
discrepancy between national laboratories, we are undertak- 
ing a thorough examination of the methods used by the other 
laboratories. The results of that work are to be reported in 
future papers. 

II. Definitions 

A. Directional-Hemispherical Reflectance 

One quantity to be measured is the spectral d/h reflec- 
tance, p (U, P, X) at a point on a plane surface. Using the 
notation described in NBS Technical Note 594-9 [10], this 
quantity can be expressed in terms of a generalized scattering 
function S as: 

p (U, P, X) = fj S (U, P; u, p, X) wda da> (1) 

where S(U, P; u, p, X) is the radiance emerging from the 
sample surface at point p in direction u due to a unit flux 
striking the sample at point P in direction U. It is assumed 
that the sample is non-fluorescent so that all of the radiation 
can be confined to a very narrow band at wavelength X. The 
integral with respect to the solid angle, to, in which the 
radiation emerges is taken over the entire hemisphere of 



directions and the integral with respect to the area a from 
which the radiance emerges is taken over the entire area from 
which the flux emerges. For the uniform, isotropic samples to 
be discussed in this paper, the d/h reflectance can be repre- 
sented by p(T, X), where T is the angle between the direction 
of incidence u and the surface normal. 

B. Directional-Hemispherical Reflectance Factor 

A typical d/h reflectometer has an entrance port which 
subtends at the sample a circularly shaped solid angle (x)' 
which in magnitude is on the order of 10~ 3 steradians. None 
of the flux falling into 0)' reaches the detector, and therefore 
the instrument does not compare d/h reflectances of the 
samples as just defined but rather compares "directional- 
hemispherical" reflectance factors [11] F(T, ,X) for collection 
over the entire solid angle except for (*)' . That is to say, the 
retroreflected flux is excluded from the measurement. To a 
very good approximation, 



F(T, A) = 



['S(U, P; 11, p, \)u • dadio 
Jf [8(P-p)/?r]u- <ia</co 



(2) 



where the prime on the integral indicates integration over the 
entire hemisphere except for a/ and 8(P — p)/vr is the 
idealized scattering function S for an ideal diffuse reflector. 
It is the quantity F(T, X) in which we are most interested 
when calibrating a reflectometer, and determining this quan- 
tity for several working standards is the object of the work 
reported in this paper. 

III. An Overview of the Method 

The method to be described for determining d/h reflec- 
tance factor is more complex than the method for measuring 
d/h reflectance originally proposed by Van den Akker. In the 
original method, it was assumed that d/h reflectance is 
independent of the direction of incidence. Since this is not 
sufficiently correct for many surfaces, additional measure- 
ments are needed in order to accurately determine the d/h 
reflectance using this approach. Also, the relationship be- 
tween d/h reflectance and the reflectance factor as measured 
by an instrument must be determined. However, for simplic- 
ity, we will continue to call the entire process the Van den 
Akker method. 

Although this method is simple in principle and is poten- 
tially one of the most accurate methods, the descriptions of it 
in the literature are usually given in terms of involved sum- 
mations of infinite series. These descriptions tend to cloud 
the simplicity of the method and to lead to a misunderstand- 
ing of it. Therefore, before describing the details of the 
measurements as performed at NBS, we will first provide a 
brief overview of the method in terms of four main steps. The 
first and last steps were used in the original Van den Akker 
method, and we have added the additional steps to take into 
account the variation of p(T, X) with T and the difference 
between p(T, X) and F(T, X). 

The primary measurement in the Van den Akker method is 
a measurement of the reflectance of the wall of an integrating 
sphere (fig 1) under the irradiation it receives in the sphere. 
We will call this reflectance the Van den Akker reflectance 
p v (k). A flux is introduced into the sphere by reflecting a 
collimated beam of radiation from the back of the sphere. 
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FIGURE 1. Sphere used for measuring Van den Akker reflectance. 



Most of the radiation which strikes a given area of the sphere 
wall has undergone many reflections. In a sphere of this type, 
the average number of times a photon will have passed into 
the sphere wall and been reflected before it is absorbed or 
escapes would be between fifty and one-hundred. Because 
the reflection from the sphere wall is very diffuse, any 
asymmetry in the way the photons are introduced is dispelled 
in one or two reflections and the distribution of radiance, L s , 
striking the wall has two characteristics properties. First, the 
irradiance is very nearly uniform over the sphere wall. This 
follows from the fact that a perfectly Lambertian (radiance 
constant with respect to direction) reflected flux would pro- 
duce a uniform irradiance over the sphere wall. Since the 
distribution of the reflected flux is very nearly Lambertian, 
the distributon over the sphere wall of incident photons which 
have been reflected more than two or three times is for all 
practical purposes uniform. This phenomenon is well-known 
and forms the basis for the many applications of integrating 
spheres in which the sphere is used as a flux averaging 
device. Another property of the spherical geometry which is 
less often noticed is that a photon which leaves the sphere 
wall at an angle 6 will next strike the sphere wall at the same 
angle. Thus, after one or two reflections, an angular distribu- 
tion of radiance L s (T), will be set up which represents, to 
within a proportionality constant, the incident and emerging 
radiance at any point on the sphere wall. Unlike the distribu- 
tion of irradiance over the wall area, the distribution of the 
radiance with respect to the angle does not become uniform 
with angle but takes on a distribution determined by the 
bidirectional reflectance properties of the wall coating. 

The remaining steps in the method relate the Van den 
Akker reflectance to the d/h reflectance factor and provide a 
method of using this reflectance data to calibrate an instru- 
ment. The four steps in the Van den Akker method can be 
outlined as follows: 

A. Determining the Van den Akker Reflectance p,(X) 

A sphere with a single port (see fig. 1) is lined with a 
highly reflecting white coating, radiation is introduced into 
the port in a collimated beam, and, from measurements made 
of the flux re-emerging from the port and of the flux reflected 



from the target area, the Van den Akker reflectance p r 2 is 
determined. A flux (j) (A.) is introduced into the sphere in a 
collimated beam which strikes the target area of reflectance 
p t (r, X), providing a nearly Lambertian source of radiation 
which emits a flux cf) p t . A fraction f of this flux passes out 
through the port as a nearly collimated beam, and the re- 
maining portion falls on the sphere wall and is reflected with 
reflectance p' v , the average reflectance of the sphere wall, 
excluding the port. This flux, the strength of which is given 
by (f>opt (1 — f)p' v , acts as a radiation source which is uni- 
formly distributed over the sphere wall. This flux is the 
source of radiant energy for the total flux </> t {X) striking the 
sphere boundary. A fraction / of Xj passes out through the 
port, and a fraction (1 — p' v ) of the remaining (1 — f) is ab- 
sorbed in the sphere wall or otherwise lost and is not returned 
to the sphere. In equilibrium, the flux lost from </>j out the 
port and into the walls is replaced by the uniformly dis- 
tributed source, from which an equilibrium flux balance 
equation can be obtained: 



<t> pti -f)p'v = <t>i\f + a -/)(i - p'v)]. 



(3) 



If the flux </> is the sample beam flux from a dual beam 
reflectometer, a signal Q J(k) is obtained with the sphere in 
the sample position: 



Q s =k[^{\ -/') + <t>oPtf(l -/")] 



(4) 



where/' is the fraction of the lambertian flux which escapes 
from the ports of the reflectometer sphere and f" is the 
corresponding fraction of the nearly collimated beam which 
emerges from the sphere as a result of the first reflection from 
the back of the Van den Akker sphere. With the target area 
from the back of the sphere placed directly on the reflectome- 
ter sample port as a sample, a signal Qi^) is obtained: 



Q t = k<f> p t (l -/'). 



(5) 



It is arranged that the angle of incidence T at which (f> 
strikes the target when Q t is measured is the same as when Q s 
is measured. Equations (3), (4), and (5) can be solved to 
obtain an expression for p' v in terms of Q t , Qs,f,f' and/". 



vher 



1 - (JQtlQMX/ (1 ~ «)) 
1-/ 



/(/" - /') Qt 

(i-f) Qs 



(6) 



(6a) 



If a were 0, this would be the usual form cited for the Van 
den Akker reflectance. The factor 1/(1 — a) takes into 
account that the first reflection of the incident beam from the 
target in the sphere emerges nearly collimated whereas the 
remaining flux emerging from the sphere is nearly Lamber- 
tian. If the reflectometer handled both fluxes in the same 
way, i.e., if/' and/" were equal, this term would be zero. 
However, it is usually the cause that/" is very much larger 
than/' and this contribution must be included. 



2 To keep the writing uncluttered, the functional dependence of each quantity is shown only when 
it is introduced and at key points in the development. 
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An additional modification must be made in the original 
Van den Akker equation to take into account the effect on p v 
of the retroreflectance of the sphere wall. This modification, 
which results in only a small adjustment p v but which greatly 
complicates the form of eq (6), is discussed in detail in 
appendix A. 

Since (1 — f)(\ — Pv) in eq (3) represents the portion of <f} t 
which is lost everywhere but out of the port, it follows that the 
reflectance p v ' is the average reflectance over the sphere area 
excluding the port. To relate p^ to a property of the sphere 
coating, it is important that the sphere coating be uniform and 
completely cover all of the sphere area except for the en- 
trance port. If it is necessary to have cracks or other gross 
imperfections in the coating, the losses in such imperfections 
must be estimated and p' v must be suitably corrected in order 
to obtain the Van den Akker reflectance p v of the sphere 
coating. 

B. Adjustment from Van den Akker Reflectance p v (\) to 
d/h Reflectance p(Y, X) 

In order that the need for the next steps be more readily 
understood, p v will be interpreted in terms of the directional- 
hemispherical reflectance, p, as defined in eq (1). 

2tt I L S (T) p(T) cos(T) sin(r) dY 
Jo 



L,(y, X) = A,(X) N (y, X)cos(y) 



(9) 



Pv 



277 l s (T) cos(r) sin(r) dr 

Jo 



(7) 



where k\ is a constant of proportionality. 

To complete the adjustment, information about the form of 
p(T, X) is obtained. A flat sample of the sphere wall coating 
is prepared in the same way as the sphere wall coating and its 
relative d/h reflectance is measured. The corrected reflec- 
tometer signal [12] N G is proportional to the d/h reflectance. 



p(T, X) = k 2 (k) A/ G (I\ A). 



(10) 



The constant k 2 can be expressed in terms of measured 
quantities by simultaneous solution of cqs (7), (9), and (10). 

, m _ pMJn.F, x) sinr dr Pv (k) 
* //v c (r,\)/v 9 (r,\)sinr</r /v G (\) ( ' 

where N G is the weighted average of N G {Y , X) with weighting 
function N g (T, X) sinT. With the value of k 2 determined, the 
d/h reflectance for any wavelength and angle of incidence can 
be calculated from the relative d/h reflectance data by using 
eq (10). However, we prefer to use an expression forp(r, X) 
which takes the form of a small adjustment of the Van den 
Akker reflectance. By substituting the expression for k 2 from 
eq (11) into eq (10) and using some algebraic manipulation, 
one can put eq (10) into the form 

p{T, \) = p„(X)|l + C(T, X.)] (12) 



where L S (T) is the radiance associated with fa at the incident 
angle T and explicit indication of the dependence on wave- 
length has been omitted for clarity. As we have indicated, p v 
is just the weighted average over all directions of incidence T 
of the d/h reflectance. The weighting function L S (T) cosl 
sinT is proportional to the irradiance of the sphere wall per 
unit angle at the angle T. If the d/h reflectance were inde- 
pendent of angle of incidence, one can see from equation (7) 
that p v and p{T) would be identical. However, in any real 
sphere coating, p(T) is not quite constant and this variation 
must be taken into account if d/h reflectance is to be deter- 
mined from p v . 

The flux, fa, striking the sphere wall is composed of 
radiation most of which has undergone many reflections. As 
has already indicated, the radiance L S (T) striking the wall in 
one location emerged-from the wall at the same angle at some 
other location. When the irradiance of the wall is uniform and 
if the generalized scattering function S of the coating were 
known, the function L s could be determined by solving the 
integral equation 



I, 



7T/2 



Pv Uy) = 277 ui) s(r, y) co S r sinr dr. (8) 

Jo 

Because the irradiance of the sphere wall is uniform, the self 
radiance distribution, L s , can be determined to within a 
constant by measuring the radiance, L s (y), emerging from the 
sphere coating through the port with an uncalibrated goniore- 
flectometer. The relative signal N g from the gonioreflectome- 
ter as a function of the angle of observation y can be related 
to L 8 (y) as 



c ( r,x) = ^-i. 

C. From d/h Reflectance p(T f X) to d/h Reflectance 
Factor F(T, X) 

If the generalized scattering function S for directions of 
incidence and reflectance both falling within the reflectome- 
ter entrance port (retroreflectance at entrance angle T) were 
equal to the average of S over all directions of viewing, then 
the reflectance p(T, X) and the reflectance factor F{T, X) will 
be equal. However, for the usual type of sphere coating, 5 in 
the retroreflective direction has been observed to be as much 
as 50 percent above the average at small angles of observa- 
tion [13]. Therefore, the bidirectional reflectance factor F(U, 
u X) must be determined over the extent of the solid angle o>' 
surrounding the incident direction U and the relationship 
between F(T, X) and p(T , X) determined from eqs (1) and (2) 
as 

p(r, X) J F(U, u)cos7 da) 

F(r,X) = - -"— (13) 



( 1 — — a)' cosy J 



where the integral in 0) is taken overo/. (Note thatF(U, u) = 
7T fS (U, u)da.) 
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D. Determining the Directional-Hemispherical Reflec- 
tance Factor of Standards and Other Samples 

The d/h reflectance factor F(T, A.) as determined in the 
preceding four steps is the d/h reflectance factor of the sphere 
wall coating. The final step is to use this information to 
measure a stable standard reflector which can then be used to 
calibrate a d/h reflectometer. One of the methods used by 
Van den Akker [9] was to compare the standard directly with 
sections of the sphere wall which were hoped to be represent- 
ative. Another method is to prepare a large number of sam- 
ples using the same materials and preparation techniques 
that were used in coating the sphere itself and comparing 
these to the standard. In either case, the d/h reflectance 
factor F C (T, A.) of the standard is calculated from the previ- 
ously determined d/h reflectance factor F(T, A) of the sphere 
wall as 



F e (T, X) = F(T, X) 



n Q C (T, X) 

n 

2 <?<(r, x) 



(14) 



where Q c is the reflectometer signal for the standard and Qi is 
the reflectometer signal for the ith one of the n samples which 
are taken to represent the sphere wall. In day-to-day use, the 
d/h reflectance factor F X (T, A) of a sample can be measured 
by comparison to the stable standard. 



F x (r, x) = F c (r, x) 



Q X (T, X) 
Qc(T, X) 



(15) 



where the reflectometer readings Q x and Q c are taken close 
together in time in order to minimize the effects of instrument 
drift. 

IV. The Measurements as Currently Made at 
NBS 

The overview in the preceding section describes in general 
the way that d/h reflectance can be determined by the Van 
den Akker method. Three determinations of the d/h reflec- 
tance of a set of standards has been made using the DRS. 
Through the experience gained in these determinations, a 
procedure for making the measurements with this instrument 
has been worked out. That procedure is described in this 
section to give a documented starting point for future im- 
provements in the method and the data from these three 
determinations provide the basis for the current NBS scale of 
6°/h reflectance factor. 

Three different materials were used as sphere coatings for 
the determinations. In the first determination, the sphere was 
coated with BaS04 powder [3, 14] from a stock which had 
been on our laboratory shelf for some time and had been 
opened on previous occasions. This was used to obtain 
experience with the method and to obtain an idea of some of 
the difficulties which might be encountered with a material 
which was slightly contaminated and which had a slightly 
lower reflectance than the best quality coatings. In the sec- 
ond determination, the sphere was coated with Halon [3, 15] 
powder, and in the third determination a fresh lot of reagent 
grade BaS0 4 powder [3, 16] was used. The advantages and 
disadvantages of each coating will be brought out in the 



discussion of the measurement errors in the description of 
each of the three main measurement steps. 

A. Determining the Van den Akker Reflectance p r (A.) 

The design of the spheres used to determine the Van den 
Akker reflectance of the sphere coating material is shown in 
figure 1. In some of his earlier work using this method, Van 
den Akker used a sphere with many removable sections in 
order to sample at many points the actual sphere coating 
being measured. However, the presence of many removable 
sections provides an uneven substructure for the sphere 
coating, creating a potential for dark rings and even cracks in 
the coating at the border of each removable section. There- 
fore, we decided to have only the target area removable and to 
rely on reproducing the coating for a representative sampling. 

The sphere coating is made by packing the powder into the 
aluminum sphere shell with an electrically driven hammer 
which has a Teflon [3] head (fig. 2). The hammer is mounted 
on a stand and each hemisphere of the sphere is mounted in 
turn on a gimbal which allows it to be rotated about a point a 
fixed distance behind the hammer. In this way, a hemisphere 
of constant interior radius R is packed uniformly in a me- 
chanical way. A plastic ring attached to the flange of each 
hemisphere as it is being packed allows the packing jig to be 
aligned the same for both hemispheres and allows the coating 
to be built up to full thickness at the edge of the hemisphere. 
In this way, when the plastic ring is removed and the sphere 
is joined, a uniform sphere coating of interior radius R results 
which has no crack at the plane where the hemispheres join. 




FIGURE 2. Apparatus used to produce the sphere coating. 



1 . Determining the Uncorrected Van den Akker Reflectance p[. 

The fraction/ introduced in eq (6) is taken as that portion 
of the area of the sphere of radius R which is occupied by the 
port of radius r, i.e., 



/=(!-(! - (r/Rf) w )/2. 



(16) 



Associated with each determinaton of/ there is an estimated 
random error, as defined in appendix B. A random error 67? 
arises from our inability to pack the sphere to the same radius 
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R from one determination of reflectance to the next and from 
the inability to measure the average/? exactly. This gives rise 
to a random error 8f which is given by: 



8/ = 



(r 2 /R 3 



2(1 - (r/R) 2 ) 



2\1I2 



8R. 



(17) 



A systematic error arises from any error in measuring r and 
from assuming that the measured r is the effective r. Both of 
these can be taken into account in terms of an uncertainty Ar 
in the effective port radius. The estimated systematic error A/ 
is given by 



V- 



(r/R 2 ) 



2(1 - (r/R) 2 ) 1 



Ar. 



(18) 



The reflectance of the sphere was measured using the DRS 
with the 6°/h (specular included) general purpose integrating 
sphere [17]. In this mode of operation, the instrument is a 
dual beam reflectometer with a capability of highly accurate 
measurements of relative 6°/h reflectance factor, i.e., the 
instrument measures a quantity Q which is proportional to the 
6°//i reflectance factor of the sample appearing in the plane of 
its sample port. With the sphere in place as a sample, a 
reflectance value which will be called Q' s is measured. As is 
described in detail in the Technical Note on the DRS [18], a 
correction must be made for the small amount of radiation 
which does not enter the sphere port but instead is reflected 
from that portion of the sphere port plate which shows in the 
reflectometer sample port or is reflected from the sphere wall 
of the reflectometer itself. In order to correct for this scattered 
radiation, a second reflectance measurement is made with the 
auxiliary sphere port plate in place without the sphere behind 
it, i.e., with most of the sample beam of the reflectometer 
passing out into the room and being lost. The reflectance 
value Q p obtained in this manner must be subtracted from Q' s 
in order to obtain a reflectance value proportional to the 
radiation being returned from inside the sphere. Most of the 
radiation emerging from the sphere passes freely back 
through the port into the reflectometer sphere. A small 
fraction of it, however, will strike the polished beveled edge 
of the auxiliary sphere port and will, thereby, suffer a slight 
loss as it returns to the reflectometer sphere. To correct for 
this loss, we measured the relative reflectance Qf of a flat 
plate of the same stainless steel from which the auxiliary 
sphere port plate was made. Since Qf in the reflectometer is 
approximately the reflectance of the sample relative to the 
sphere wall which in turn has a reflectance nearly 1, Qf was 
approximately the reflectance of the stainless steel plate for 6° 
incidence. The radiation coming from the sphere which 
strikes the beveled edge strikes it at near grazing incidence. 
Since the reflectance from a metal surface near grazing 
incidence is higher than the near-normal reflectance, we 
chose to represent the reflectance from the beveled edge by (1 
+ 2Qf)l3- This estimate is based on the general shape of the 
reflectance curves for metals [19] and on the port dimensions 
(fig. 3) [20]. Under this assumption, if the projected solid 
angle subtended by the beveled edge of the port averaged 
over the port area is 0> e , an expression for the corrected 
sphere reflectance Q s can be written as 



<?, = (Q's - Qp)l\ i 



w e 2(1 - <?. 



*]■ 



(19) 



POLISHED BEVELED PORT EDGE 
(STAINLESS STEEL) 




Figure 3. Cross section detail qf polished beveled port edge. 



Associated with this determination of Q s , there is a random 
error 8Q 8 arising from the random errors in measuring Qs, 
Q p , and Qf. Because of the smallness of G) e , the denomina- 
tor in eq (19) is approximately 1 and the small random error 
in Qf has little effect on the total random error. Therefore 



8Q S ^(8Q' ! 2 + 8Q 2 y 



(20) 



Also associated with the determination of Q s is a systematic 
error AQ S which arises from the uncertainty in Q) e and from 
the approximation used in estimating the effective reflectance 
of the beveled edge of the port. The systematic uncertainty in 
the beveled edge correction is estimated to be one-half as 
large as the correction itself, so that 

2(1 - Qf) ( /o>, \ 2 \ 1/2 

The reflectance of the removable target was measured 
using the same instrument. In this case a value Q\ was 
obtained which also included a small contribution due to 
stray radiation from the sample beam which is reflected from 
the reflectometer sphere wall and from that portion of the 
target which is outside the region occupied by the auxiliary 
sphere port. A correction for this contribution was obtained 
by measuring Q e with a sample in the port made of the same 
material as the target but which has a hole in it the size of the 
port in the auxiliary sphere. Note that in this measurement, 
the beveled edge is not present, since the stainless steel port 
structure depicted in figure 3 is part of the auxiliary sphere. 
The corrected relative reflectance of the target Q t is obtained 
as 



Q, = Q'<- Qe- 



(22) 



An expression for the random error associated with this 
determination of Qt is determined in a straightforward fashion 
and resembles eq (20). 



SQt = W, 2 + SQW 



(23) 



The values of/, Q s , and Q t obtained above were used in eq 
(6) (more specifically, eq (A7)) to calculate values of p^X), 
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the Van den Akker reflectance of the sphere wall in its 
entirety. The random error in p v is obtained from eq (6) as 



8p; = [(d -Qt/Q,)l(i -/)W 

+ (if/Q s )/a -f)?m 

+ WQtlQDia -f)f apff*. 

The systematic error in p' v is given similarly by: 

Api=[((l -QJQ.)I(\ ~J) 2 ) 2 Af 



(24) 



(25) 



where AQ S is only that portion of the systematic uncertainty 
in Q s given by eq (21). Note that the contribution to the 
systematic error in Q s and Q t due to systematic errors in the 
DRS are omitted because these errors are very small [21] and 
because they tend to cancel when the ratio QJQs is taken in 
determining p v from eq (6). Therefore there is no AQ t term in 
eq (25). 



2. From Uncorrected Van den Akker Reflectance p' v to Corrected Van 
den Akker Reflectance p v 

The quantity p' v as determined in the preceding section is 
the average Van den Akker reflectance of the sphere wall. In 
order to obtain the Van den Akker reflectance of a thick 
coating of the sphere wall material, it is necessary to correct 
q' v for the effect of the crack around the removable target 
section and for translucency in the sphere wall coating. 

a. The Gap Around the Target 

Since the target portion of the sphere is removable, there is 
a narrow gap in the sphere wall around the target. The gap 
acts as a light pipe between two aluminum surfaces, so that 
the effective reflectance of the gap is essentially zero. There- 
fore, in order to relate the measured Van den Akker reflec- 
tance p' v to the Van den Akker reflectance p v of the wall 
coating, the loss in the gap must be accounted for: 



Pv 



Pv 



k t [l-2r t W/(4R 2 -r*)] 



(26) 



where r t is the radius of the target, W is the width of the gap, r 
is the radius of the entrance port, and ^ is a transluency 
correction factor as determined in section IV. A.2.b. below. 
The systematic error associated with the crack correction is 
estimated to be one-half the magnitude of the correction. The 
total systematic error in p v is given by 



Ap r = p v 



Ap» 
P 



tr+(t) 



r t W 



4/? 2 - r 2 



(27) 



are depicted in figure 4. The coating in the sphere is, on the 
average, about 4 mm thick, so that the reflectance of a 
barium sulfate sphere wall is equal to the reflectance in the 
limit of a very thick wall, p^ [22]. Therefore, for these 
coatings k t — 1. However, in the case of Halon, the wall is 
not thick enough and a correction was made based on the data 
in figure 4. In this correction, the thickness of the coating at 
any point was taken to be the distance from the surface to the 
aluminum substrate at that point, and a translucence correc- 
tion factor k t was obtained by a calculation having the follow- 
ing form: 



k t = 



Pio 
JdA 



(28) 



where p x is the relative reflectance for the coating thickness 
at a given point, 

Pio is the relative reflectance for a 10 mm thick coating 
(assumed equal to pj, 

dA is an element of the sphere wall area, 
and the integral is taken over the entire area of the sphere 
wall. 
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FIGURE 4. Reflectance of wall coatings as function of thickness. 

(Normalized to give a reflectance of approximately 1 for thick samples.) 



b. Translucency of the Wall Coating 



The relative reflectance as a function of thickness was 
measured at each of three wavelengths on samples of each 
type of coating material. The results of these measurements 



It is estimated that the coating thickness measurements 
might have been in error by as much as V2 mm for the 4 mm 
thickness. Therefore the calculation of eq (28) was redone 
using 3.5 and 4.5 mm for the average coating thickness over 
the main sphere wall in order to obtain an estimate of the 
random uncertainty dk t of this correction. An expression for 
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the random error in p v as given by eq (26) is 



&Pv = Pv 



Pv) \kt) . 



(29) 



Since using the distance from the coating surface to the 
aluminum substrate as its thickness is only a rough approxi- 
mation in the regions of the port and of the target boundary, a 
systematic error will be introduced by that approximation. An 
estimation of this uncertainty was obtained by calculating k t 
= P4/P10, and estimating the systematic error as (k' t — k t )/2. 
We assumed the translucency correction to be approximately 
independent of wavelength. The data in figure 4, for 1 mm 
thickness indicates that this may not be exactly correct. 
Therefore, there may be a wavelength dependent systematic 
error of approximately ±0.0005 in k t , with the positive error 
associated with short wavelengths and the negative error with 
long wavelengths. 

The Halon coatings used in our apparatus had a density of 
0.8 g cm -3 . Grum [23] refers to a coating of Halon which is 
opaque at 2 mm thickness. However, from the description of 
the pressure under which the coating was formed, his coat- 
ings are probably much denser than ours. 

3. Sample Curvature and Sphere Size 

Questions have been raised concerning the effect of the 
size of the sphere upon the Van den Akker reflectance. 
Possible errors may arise due to the inability of the measuring 
instrument to properly compare the reflectance of curved and 
flat surfaces, either due to the difference in the relative areas 
of the port and the remainder of the sphere wall or due to a 
change in the reflectance of the wall coating with curvature. 
This question has been addressed experimentally in two 
different ways. The relative 6°//i reflectance of curved and 
flat samples was measured directly. The average relative 
reflectance at 550 nm of four concave BaS0 4 [14] samples 
from the first determination was 0.9826 ± 0.0027 and the 
average reflectance of two flat samples of the same material 
was found to be 0.9834, and the uncertainty is probably of 
the same order of magnitude. The difference of 0.0008 is 
probably not significant. 

In the second test, the Van den Akker reflectance was 
determined using two different sized spheres, one with a 95 
mm interior radius and the other with a 146 mm interior 
radius. Both spheres were coated with the type of BaS04 [6] 
used in the third determination. The Van den Akker reflec- 
tances for these two determinations are given in table III. It 
can be seen that there is no significant difference between the 
results from the two spheres. 

4. Results and Error Analysis 

The Van den Akker reflectances p v determined for each 
coating at twenty-five nanometer intervals over the wave- 
length range 400 nm to 750 nm are given in tables IV, V and 
VI, and are pictured in figure 5. A listing of typical values 
used in the calculations for the Van den Akker reflectance is 
shown in table VII, and representative calculations of the 
random and systematic uncertainties in p v are presented in 
tables VIII and IX respectively. The only known wavelength 
dependency of these uncertainties is in the random errors. 



Table III 

Comparison of the Van den Akker reflectance of p v of the walls of two 
different size spheres coated with BaS0 4 | 16]. 



Wavelength 


Sphere 


radius H 


DifT.. ...,,.... ;,, n 


(nm) 


14.6 cm 


9.5 cm 




400 
425 
450 
475 


0.9744 
.9771 
.9789 
.9802 


0.9744 
.9769 
.9788 
.9801 




0.0002 
.0001 
.0001 


500 
525 
550 
575 


.9813 
.982.3 
.9829* 
.9834 


.9812 
.9821 
.9828 
.9832 


.0001 
.0002 
.0001 
.0002 


600 
625 
650 

675 


.9836 

.9838 
.9838 
.9837 


.9835 
.9836 
.9836 
.9836 


.0001 
.0002 
.()(X)2 
.0001 


700 
725 
750 


.9836 
.9835 
.9832 


.9835 
.9833 
.9831 


.0001 
.0002 
.0001 



For this value, €p v = 2.1 X 10 4 . 

Table IV 
First Determination BaS() 4 |14] 



Wavelength 
(nm) 



Van den Akker 
Reflectance p v 



P(6°, A) 



F(6\ A) 



400 


0.9661 


0.9606 


0.9599 


425 


.9688 


.9638 


.9031 


450 


.9721 


.9675 


.9668 


475 


.9744 


.9703 


.9696 


500 


.9767 


.9728 


.9721 


525 


.9788 


.9752 


.9745 


550 


.9806f 


.9774* 


.9767@ 


575 


.9823 


.9793 


.9786 


600 


.9835 


.9807 


.9800 


625 


.9845 


.9818 


.9811 


650 


.9852 


.9826 


.9819 


675 


.98.58 


.9834 


.9827 


700 


.9864 


.9840 


.9833 


725 


.9869 


.9846 


.9839 


750 


.9873 


.9850 


.9843 



t ep v is 2.3 X 10~ 4 (8p v = 2.1 X 10~ 4 and Ap r = 1.1 X 10 -4 ). 
* ep is 9.8 X 10" 4 (8p = 9.5 X 10~ 4 and ^p = 2.3 X 10" 4 ). 
@e8 is 1.01 X 10~ 3 {8F = 9.5 X 10~ 4 and AF = 3.3 X 10~ 4 ). 

Table V 

Second Determination Halon [15] 



Wavelength 
(nm) 



Van den Akker 
Reflectance p v 



P(6°, A) 



F(6°, A) 



400 


0.9961 


0.9930 


0.9927 


425 


.9961 


.9932 


.9929 


450 


.9961 


.9933 


.9930 


475 


.9960 


.9934 


.9931 


500 


.9959 


.9935 


.9932 


525 


.9958 


.9935 


.9932 


550 


.99571 


.9935* 


.9932@ 


575 


.9956 


.9935 


.9932 


600 


.9955 


.9935 


.9932 


625 


.9954 


.9935 


.9932 


650 


.9952 


.9934 


.9931 


675 


.9951 


.9934 


.9931 


700 


.9950 


.9933 


.9930 


725 


.9949 


.9932 


.9929 


750 


.9948 


.9931 


.9928 



t ep v is9.0 X H)- 4 (Sp r = 7.3 X 10 ~ 4 and Ap v = 5.1 X 10 ~ 4 ). 
* epis9.6 X 10 -4 (6p= 7.9 X 10 ~ 4 and Ap = 5.5 X 10" 4 ). 
@eF is 9.7 X 10" 4 (8F = 7.9 X 10~ 4 and AF = 5.6 X 10~ 4 ). 
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Table VI 



Third Determination BaS0 4 [16] 



1.00 



Wavelength 
(nm) 


Van den Akker 
Reflectance p v 


P(6°, A) 


F(6°, A) 


.99 


400 


0.9744 


0.9695 


0.9688 


<C 


425 


.9771 


.9732 


.9725 


UJ 


450 


.9789 


.9757 


.9750 


u 

z 


475 


.9802 


.9775 


.9768 


g .98 


500 


.9813 


.9791 


.9784 


u 


525 


.9823 


.9804 


.9797 


_j 


550 


.9829t 


.9813* 


.9806@ 


UJ 


575 


.9834 


.9819 


.9812 


QC 


600 


.9836 


.9822 


.9815 


QC 


625 


.9838 


.9824 


.9817 


111 


650 


.9838 


.9824 


.9817 


* ( 


675 


.9837 


.9822 


.9815 


< 


700 


.9836 


.9819 


.9812 


2 .96 


725 


.9835 


.9815 


.9808 


Q 


750 


.9832 


.9808 


.9801 


Z 



l ' ' H 

X X X X X X X X X X X X X x * 



o o °° 



t ep v is 2.3 X 10~ 4 (8p v = 2.1 X 10~ 4 and kp v = 1.1 X 10~ 4 ). 
* ep v is 2.1 X 10~ 3 (dp = 2.1 X 10~ 3 and Up = 2.3 X 10~ 4 ). 
@eF is 2.1 X HT 3 (8F = 2.1 X 10~ 3 and AF = 3.3 X 10~ 4 ). 



The smallest random error is in the central region of the 
spectrum, where the product of the source intensity and the 
receiver sensitivity is maximum. At the short wavelength end 
of the spectrum, there is an increase in uncertainty in Q t and 
Q s associated with a decrease in source intensity, while at the 
long wavelength end of the spectrum there is an increase in 
uncertainty associated with a decrease in receiver sensitivity. 
However, because the random error associated with/ domi- 
nates the overall uncertainty, we will cite only error figures at 
550 nm for p v as representing the entire spectrum. 

The total uncertainty in determining the Van den Akker 
reflectance is obtained by adding the random and systematic 
uncertainties in quadrature. The total uncertainty in the Van 
den Akker reflectance is 0.0002 for the BaS04 coating. It can 
be seen that this uncertainty is very small, confirming the 
results of the error analysis by Goebel, et al., [2]. The larger 
uncertainty for the Halon coating is due to the transluscence 
of the coating and is introduced in the extrapolation to the 
reflectance of a thick enough layer rather than being caused 
by an uncertainty in the measurement. This additional uncer- 
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Figure 5. Van den Akker reflectance of three coatings. 

(The height of the symbols X represent the total uncertainty. Height of the symbols and A 
represents approximately 4 times the total uncertainty.) The associated solid curves show p(6°, A) for 
each coating. 

tainty would not exist if the Halon coating were twice as 
thick. 

B. Adjustment from Van den Akker Reflectance p r (X) to 
6°/n Reflectance p(6°, k) 

The second main step is the transition from p^A.) to p(6°, 
A). As was pointed out in section III, this step is an intrinsic 
part of the determination which is needed because the sphere 
coating is a real material and not an ideal Lambertian reflec- 
tor. For this reason it is not proper to regard this step as a 
correction. However, for the type of sphere coatings being 
examined, the departure from Lambertian reflection is small 
and the difference between p v and p is small. For this reason, 
we will refer to the transition step from p v to p as an 
adjustment in order to emphasize the small size of the change 
with respect to the quantity being changed. 



Table VII 
Typical Values Used in Calculating the Van den Akker Reflectance p v (A. — 550 nm) for BaS0 4 [14]. 



Symbol 



Reference Equation 



Value 



Comments 



(a) 


r 


(b) 


R 


(c) 


f, 


(d) 


Qs 


(e) 


Qp 


d) 


Qr 


(g) 


Q) e 


(h) 


Qs 


(i) 


Qt 


(!) 


Qe 


(k) 


Qt 


(1) 


Pv 


(m) 


r t 


(n) 


w 



(p) Pv 



(16) 


12.70 mm 


(16) 


146 mm 


(16) 


1.90 X 10" 3 


(19) 


8.731 X 10~ 2 


(19) 


1.24 X 10~ 3 


(19) 


7.0 X 10" 1 


(19) 


3.78 X 10~ 2 




radians 


(19) 


8.630 X 10~ 2 


(21) 


9.832 X 10" 1 


(21) 


1.36 X 10" 4 


(21) 


9.831 X 10- 1 


(6) and (A7) 


9.805 X 10" 1 


(26) 


25.4 mm 


(26) 


0.17 mm 


(26) 


1 




(0.9974) 


(26) 


9.806 x 10 _1 



measured with inside caliper and micrometer, 
measured with steel rule, 
calculated using r and R above, 
measured reflectance of auxiliary sphere, 
measured with stainless port open to dark room, 
measured reflectance of stainless steel port plate, 
calculated from dimensions of port structure [17]. 

calculated using Q' s , Q p , Qf, and (O e above. 

measured reflectance of target. 

measured with reflectometer port open into dark room. 

calculated using Q' t and Q e above. 

calculated using/, Q s and Q t from lines (c), (h) and (k) 

above respectively, 
measured with micrometer caliper, 
measured with traveling microscope. 
BaS0 4 

(Halon determination only), 
calculated using p' v , r t , W, and kt above. 
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Symbol 



Table VIII 
Calculation of the Random Error in the Van den Akker Reflectance 



Reference Equation 



Value 



Comments 



(a) 8R 



(b) 


¥ 


(c) 


SQ's 


(d) 


8Q P 


(e) 


8Qs 


(f) 


6<?; 


(g) 


8Qe 


(h) 


8Q ( 


(i) 


Sp'v 


(J) 


8k t 


(k) 


8p v 



(17) 

(17) 
(20) 
(20) 
(20) 
(23) 
(23) 
(23) 
(25) 

(27) 

(27) 



2 x 1(T 5 

7 x 10~ 5 

7 x 1(T 5 

1 x 10~ 4 

2 x 10~ 4 
7 X 1(T 5 

2.1 X 10~ 4 

2.1 x 1Q~ 4 

(7 x 1(T 4 )* 

2.1 x 10- 4 

(7.3 x 10~ 4 )* 



Uncertainty in packing sphere coating to constant ra- 
dius. 

Calculated from 8R above. 

Uncertainty from 3 measurements. 

Roundoff uncertainty. 

Calculated from 8Q' S and 8Q P as given above. 

See comment on (d) above. 

See comment on (d) above. 

Calculated from 8Q' t and 8Q e above. 

Calculated from 8f, 8Q S , and AQ t from (c), (g) and (j) 
above. 

Corresponds to an uncertainty in coating thickness of 
0.5 mm. 

Calculated from 8p' v and 
8k t above. 



* Halon coating only. 



Symbol 



Table IX 
Calculation of the Systematic Error in the Van den Akker Reflectance 



Applicable Equation 



Value 



Comments 



(a) Ar 



(b) 


A/ 


(c) 


Ao) e 


(d) 


AQ. 


(e) 


Ap; 


(0 


Ak t 


(g) 


r t W 




(R 2 + r 2 ) 


(h) 


Ap„ 



(18) 

(18) 
(21) 
(21) 

(25) 

(27) 

(27) 
(27) 



0.03 mm 

9 X 1(T 6 

2 X 10" 3 rad 

1.1 X 10~ 4 

9.7 X 10" 5 

(5 X 1(T 4 )* 

5 X 10- 5 



1.1 X 1(T 4 
(5.1 X 10" 4 )* 



Uncertainty in measuring port ±.01 mm. Another 
contribution is uncertainty in amount scattered back 
from port lip, which is estimated. 

Calculated from Ar above. 



Calculated from Aoi e above and from Qf, (D e , and Q 8 
from entries (f), (g) and (h) in table VII. 

Calculated from A/ and AQ S from (b) and (d) above. 

Wavelength dependent error ±.0005 estimated from 
fig. 4; uncertainty due to tapered edges of coating 
±0.0001. 

Calculated using/?, r, r t , W from (a), (b), (m), and (n) 
in table VII. 



Calculated using Ap ' v , Ak t , and the gap uncertainty as 
given in (e), (f) and (g) above respectively. 



Halon coating only. 



1. Determining the Relative Self-Radiance N g {y)/cos y 



For each type of sphere coating, the relative self-radiance 
N g (y)/cos y was determined as a function of angle of emer- 
gence. A 75 mm radius sphere was coated using the same 
procedure as was used when coating the spheres for the Van 
den Akker reflectance determinations. This sphere was irra- 
diated through its single port and the flux emitted from the 
port past a knife-edge was observed at a number of angles. 
The geometry of this measurement is illustrated in figure 6. 



The instrument used was the NBS-Gaertner goniophotometer 
illustrated in figure 7. The source aperture permitted colli- 
mation to within 0.25° of the optical axis, and the receiver 
was collimated to within 0.6° of the optic axis. All three 
coatings were measured with a Kx"iH urnmant C [24] spectral 
weighting, the centroid of which falls at approximately 550 
nm. In addition, the coatings from the second and third 
determinations were measured using a 550 nm low pass 
interference filter, which resulted in a centroid of spectral 
weighting at 450 nm, and using a 600 nm high pass filter 
which resulted in a centroid spectral weighting at approxi- 
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COLLIMATED 
RECEIVER 




COLLIMATED 
BEAM 



COLLIMATED 
SOURCE 



FIGURE 6. Arrangement of knife-edge apparatus for measuring N (y). 



mately 700 nm. No distinction between the measurements 
under different spectral weightings were observed, indicating 
that Lj^y) depends principally on the geometrical nature of 
the surface. The relative magnitude of the observed flux 
divided by the cosine of the observation angle is shown in 
figure 8. Also shown in figure 8 are measurements made by 
scanning, with a telescopic detector, the interior of a 20 cm 
radius sphere coated with Halon. These latter measurements, 
believed to be much more accurate (±0.001) than the knife- 
edge measurements, extend only to 45°. Because of the 
method used to produce the coating by hammering, there is 
reason to expect an increase in N g (y)/cos (y) as y approaches 
90° because of a slight glossiness in the surface. However, it 
appears that this effect is not nearly as great as the measure- 
ments using the knife-edge indicate. Two possible sources for 
the difficulties at large y are reflections from the corner of the 
knife-edge and scattered light from the receiver optics. At 
large y, the flux to be measured is small, so that small 
amounts of stray radiation can cause large errors. Fortu- 
nately, this uncertainty in NJy)/cos y results in only a small 
uncertainty in determining C(T, A.). A discussion of our 
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FIGURE 8. Normalized self-radiance at 550 nrn for each of three coating 
materials. 



choice of function to represent NJ,y) will be deferred to part 
IV. B. 3 in which C(T, k) is calculated. 

2. Determining the Relative Directional-Hemispherical Reflectance 

N G (r) 

The relative directional-hemispherical reflectance of a 
sample of each of the three coatings was measured with the 



PH0T0MUTIPLIER 



PERPENDICULAR 




RECEIVER 
APERTURE 



'stray FLUX 

SHIELDS 
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RECTANGULAR STOP 
8 SHUTTER 



RECEIVER DIAPHRAGM 



FIGURE 7. NBS-Gaertner goniophotometer. 
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DRS using the special measurement accessory apparatus for 
this purpose [12]. A diagram of this measurement accessory 
is reproduced from that reference in figure 9. The reflectance 
values were measured for radiation polarized with the electric 
vector in the plane of incidence and perpendicular to the 
plane of incidence. The average of the relative reflectance for 
the two polarizations was determined as a function of angle of 
incidence for each material for each of three wavelengths, 
450 nm, 550 nm, and 750 nm. The data were fitted with a 
quadratic function, using a least-squares fitting program, to 
obtain the /V^T, X) shown in figure 10. These curves were 
used to determine the adjustment function C(T, A.) of eq (12). 
The general trend is for the reflectance to increase as the 
angle of incidence approaches grazing. The amount of the 
increase is greater for the materials with lower reflectance. 



1.015 



SAMPLE 
SUSPENSION ROD 



ENTRANCE PORT 
APERTURE 



BEAM SIZE: 

10 MM SQUARE 

FOR 10-NM 

BAND-PASS 

(HEIGHT ADJUSTABLE) 




SAMPLE TURNTABLE 
(STEPPING MOTOR CONTROLLED) 
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PLEXIGLAS ROD 



PHOTOMULTIPLIER 
DETECTOR 



SHUTTER 
PHOTO -CATHODE 



DIRECTIONAL - HEMISPHERICAL 

REFLECTANCE 

(45 CM DIA. SPHERE) 

FIGURE 9. Integrating sphere for measuring relative directional-hemi- 
spherical reflectance as a function of angle of incidence. 



3. Determining C(T, k) 

For each type of coating, the function C was evaluated 
using in eq (13) the NJy) and Nr^Y, X) data from parts 1. and 
2. above for each of the three wavelengths 450, 550, and 750 
nm. Actually three sets of C were calculated for each of three 
interpretations of the highly uncertain data for N J,y) in figure 
8 above in order to evaluate the effect of that uncertainty. We 
wish to distinguish between the data taken from measure- 
ments of the radiance from the sphere wall and the interpreta- 
tions of this data which are used as incident radiance values 
in calculating C(T, X). We will do this by using a lower case 
y in Nj,y)/cos y to indicate the data obtained directly from 
the measurements and by using an uppercase T in N^T)/ cos 
r to indicate the processed data used as the relative radiance 
distribution incident on the sphere wall when we calculate 
C(T, X). The first set of C(T, A.) was calculated assuming that 



l« f RUN 

(B0SO4 ) 




30 60 90 
ANGLE OF INCIDENCE T (DEG.) 

FIGURE 10. Relative directional- hemispherical reflectance N g (V, K) as a 
function of angle of incidence T for each of three coatings at three wave- 
lengths. 



/Vg(y)/cos y is independent of y. This provides a lower limit 
on C. The second was calculated by using 1 for/V^y)/cos y 
between and 45°. To obtain the /V^y)/cos y above 45° in 
this approximation, a straight line was drawn through the 
knife-edge experiment data from to 45° and the amount by 
which the experimental data for greater than 45° fell above 
this line was added to I to obtain /V ^y). (See fig. 1 1 for a 
diagram illustrating this procedure for the Halon data in fig. 
8.) The set of /V^y) obtained in this way were consistent with 
the more accurate data from the telescopic detector scan and 
the C(T, X) calculated with these /V J^y) are used as the 
accepted values. Finally, as an upper limit, Cs are calcu- 
lated using the unmodified knife-edge experiment data given 
in figure 8. The functions C(T, X) corresponding to the 
accepted values, are plotted against T for each of the three 
materials at each of three wavelengths in figure 12. (See 
appendix C for details of the calculations leading to this 
figure.) 

Since the goal of this section is to determine p(6°, A.) at 25 
nm intervals over the wavelength range 400 nm ^ a < 750 
nm, we need values of C(6°, X) at these wavelengths. Since 
the large amount of data required would make impractical 
determining all of these values of C(6°, X) in the way de- 
scribed above, we choose to determine C(6°, X) from the data 
at the three wavelengths 450, 550, and 750 nm by interpola- 
tion and extrapolation. This determination is based on the 
definition of C(6°, X) which comes from eq (12) 



C(6°, X) = 



/V G (6°, X) 



1. 



(30) 
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FIGURE 11. An illustration of the three assumptions used in interpreting 
the N g {y) data from the knife-edge determinations. 
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550nm 
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Figure 12. C(T, A.) as a function of T at three wavelengths for each of 
three coatings. 



Making use of this, we can write a correction factor C'(6°, X) 
which is a rough approximation to C(6°, X) as: 



„ „ x N G (6°, X) 



(32) 



Since C and C' are approximately equal, their ratio is a 
smoothly varying function of wavelength. Therefore, we de- 
termined J(k) as 



/(X) = d + d t k + d 2 k 2 



such that 



JW = 



C(6°, X) 
C'(6°, X) 



(33) 



(34) 



at the wavelengths 450 nm, 550 nm and 750 nm. We then 
determined C(6°, X) at other wavelengths as 



C(6°, X) = 7(X)C'(6°, X) 



(35) 



where C'(6°, X) is calculated from expression (32) using 
values of N G (6°, X) and /V G (45°, X) measured with the relative 
d/h reflectance instrument [12]. The values of C(6°, X) so 
determined are plotted in fig. 13. 



-.001 



-.002- 



-.003 -> 



UJ o^ 
O 



-.004 



o o 

u 
a: 



-.005 



1 1 






/ 


N 


. — 


/ ^ --" 




\ 


/ ^^ 




\ 


' y 






A . 






^7 y^ 






_ / ! / 






' ' / 






1 / 

- l / 
i / 










1 / 










- 






/ 2"d RUN (HALON) — 







/ 3rd RUN(BaS0 4 ) 

1 1 







400 



500 600 

WAVELENGTH (nm) 



700 750 



FIGURE 13. C(6°, A) as a function of wavelength for each of three 
coatings. 



4. Results and Error Analysis 



Since the weighing to determine N o{k) is heaviest at 45 , it 
follows that 



/V G (X) - /V G (45°, X). 



(31) 



The three sets of values forC(6°, X) determined above were 
used in eq (12) to calculate p(6°, X) for the three wall 
materials. The results of this calculation are shown in the 
third columns of tables IV, V and VI and by the line curves in 

fig. 5. 
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From eq (12) one can determine a propagation of error 
formula for the random error 8p(6°, A) as: 

8p(6°, A) = [Sp„(A) 2 + 8N G (6°, kf + 8N c (k) 2 ] 112 . (36) 

In this form p r , N G and /V^are each assumed to be 1, and 
errors common to N G and N G are ignored, since these will 
effectively cancel when the ratio is taken. Therefore the only 
errors to be included in eq (36) are the independent errors in 
the various quantities. The principal independent uncertainty 
in N G is that due to the uncertainty in the self-radiance 
measurements as shown in appendix C. The principal random 
uncertainty in N G is due to the noise in the d/h data. Since 
the value of /V G (6 , A.) was obtained through a rather indirect 
procedure involving curve fitting, it is difficult to establish a 
good theoretical basis for the error estimate. Therefore, we 
sought a reproducible way of estimating the error which 
depended as little as possible on arbitrary judgment. Since 
most of the curves involved fitting through ten to fifteen 
points with a quadratic function with two independent param- 
eters, there is almost no likelihood that a point taken from the 
curve will depart from the most probable value by more than 
the root mean square deviation of the individual points with 
respect to the curve. Therefore, we will use this deviation as 
the estimated random uncertainty. In determining N G , two 
such fittings are involved, one to the original data as a 
function of angle and one to the C(A) data used in the 
interpolation. A summary of the random error analysis for the 
p v {k) to p(6°, A.) adjustment is given in table X. 

The only source of systematic error which we have identi- 
fied for this adjustment step is an uncertainty in the measure- 
ment of N G due to the correction which is made for the 
entrance port of the reflectometer. This results in an uncer- 
tainty in C(T, A.) of 2 X 10 -4 independent of wavelength. « 

The total uncertainty in p(6°, A) for a given determination 
is obtained by adding 5p(6°, A) and Ap(6°, k) in quadrature. 



The additional uncertainty introduced in making the adjust- 
ment from p v (A) to p(6°, A) is of the same order of magnitude 
as the uncertainty in p v (A) itself. Particularly noteworthy is 
the large increase in uncertainty in the third determination. 
The barium sulfate used in this determination had a rather 
"sticky" consistency which made it difficult to pack it into a 
uniform smooth surface, and this may have had an effect 
upon the noise in the N G (T, A) data from this surface. 

C. Adjustment from d/h Reflectance p(6°, A) to d/h 
Reflectance Factor F(6°, A) 

In order to determine the bidirectional reflectance factor in 
the retroreflective direction, apparatus was set up as custom- 
arily used for measuring the coefficient of luminous intensity 
of retroreflective sheeting (See fig. 14). Since the retroreflect- 
ance of the sphere coatings is essentially spectrally nonselec- 
tive over the visible wavelength range [13] we measured the 
luminous reflectance factor for CIE illuminant A as being 
representative of the entire spectrum to within the uncertainty 
of the retroreflectance measurements. 



SAMPLE SURFACE 




PROJECTOR 



Figure 14. Diagram of apparatus used in retroreflectance measure- 
ments. 



Table X 

Calculation of the errors in the p„(X) to p(6°, X) adjustment (See equation (36) and the accompanying discussion). 



Symbol 



BaS0 4 [14] 



Halon [15] 



BaS0 4 [16] 



Comments 



(a) 


8p v 


(b) 


sa/ g (X) 


(c) 


oW fl (6°, X) 


(d) 


syv,(6°, \) 


(e) 


SA/ fl (6°, X) 


(0 


Sp(6°, X) 


(g) 


&Pv 


(h) 


AC 


(i) 


Ap(6°, X) 



2.1 x 1(T 4 


7.3 X 


10" 4 


2.1 X 


10" 4 


See Table VIII. 


0.5 X 10" 4 


0.5 x 


io~ 4 


0.5 x 


10" 4 


See Appendix C. 


7 x 10" 4 


3 X 


10" 4 


20 X 


10" 4 


Component due to fitting angular data 
for/V G (r, X). 


6 x 10" 4 


3 x 


10" 4 


5 X 


10" 4 


Component due to fitting wavelength 
data to C'(X). 


9.2 x 1(T 4 


4.3 x 


10" 4 


21 x 


10" 4 


Quadrature combination of (c) and (d) 
above. 


9.5 X 1(T 4 


7.9 x 


10" 4 


21 x 


10" 4 


Quadrature combination of (a), (b) and 
(c) above. 


1.1 x l(T 4 


5.1 X 


10~ 4 


1.1 x 


10" 4 


See Table IX. 


2 x 10" 4 


2 x 


10~ 4 


2 x 


10" 4 


Associated with d-h reflectometer port 
correction. 


2.3 x in- 4 


5.5 X 


10" 4 


2.3 X 


10~ 4 


Quadrature combination of (g) and (h) 



above. 
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To make the measurements, the receiver was placed in the 
sample position and a signal Nr proportional to the normal 
illuminance I r on the sample was measured. 



N* 



^rIr- 



(37) 



The receiver was then moved to a position a distance d away 
from the sample and, on the same scale of measurement and 
using the same receiver aperture, signals A^/3), which is in 
the same way proportional to the illuminance I(J3) on the 
receiver due to the radiation reflected from the sample, was 
measured with the entire sample area A in view. From the 
basic definition for reflectance factor, it follows that the 
reflectance factor F r (/3) can be determined from 



F r (P) = 



NAfi) 



n r d 2 cos (D cos (r + /3)' 



(38) 



The value ofF 6 o(/3) so determined as a function of observation 
angle fi is given for BaS0 4 and for Halon in fig 15. (See 
appendix A, sec. A. 2 for further discussion of these data and 
the evaluation of the integrals.) The retroreflectance factor 
F 6 °(/3) was found to be essentially independent of the angle 
which the plane of observation makes with respect to the 
plane of incidence. Therefore, the integral in eq (13) be- 
comes 

J*"' F(U, u) cos y d <o = 2n cos6° J r ld 'F 6 lp)f3d p (39) 

and O)' in (13) becomes 

«u' = 2TTj r ld 'l3dl3 (40) 

where r is the radius of the reflectometer entrance port and d' 
is the distance in the reflectometer sphere from the sample to 
the plane of the entrance port. The p(6°, X) data were ad- 
justed to F(6°, X) using eq (13). In general, the difference 
p(6°, X) - F(6°, X) is 0.0007 for BaS0 4 and 0.0003 for 
Halon. The estimated systematic uncertainty in this adjust- 
ment is approximately one-third of its value in each case. 
These uncertainties, added in quadrature with the uncer- 
tainties in p(6°, X) yield the total systematic uncertainty in 
F(6°, X). The values of F(6°, X) for the three coatings are 
given in the third columns of tables IV, V, and VI. 

D. Determining the d/h Reflectance Factor of Working 
Standards 

The d/h reflectance factor F (6°, X) determined in IV. C. 
above is the average for the wall coating used in the Van den 
Akker sphere. Since this coating is neither permanent nor 
readily accessible, the reflectance factor F c (6°, X) is deter- 
mined for working standards made of glass, ceramic, or some 
other permanent material. The value of F c is assigned as 
indicated in eq (14) and the accompanying text. Specifically, 
we will denote by F c j(6°, X) the reflectance factor for stan- 
dard c as obtained from theyth determination. Equation (14) 
takes the form: 



ft ' 



(41) 



where 

Fj =Fj(6°, X) is the reflectance factor for the yth wall 
coating, 



/J*f 



6° ANGLE OF INCIDENCE 

1st BaS04 
X 2nd Halon 




OBSERVATION ANGLE 

Figure 15. Retroreflective bidirectional reflectance factor for 6° inci- 
dence as a function of observation angle. 



Q c is the measured relative reflectance of working stan- 
dard c, and 

Qj is the average of the measured relative reflectance for 
several representative samples of they'th wall coating. 

The relative reflectances Q c and Qj were found by experiment 
to be independent of the polarization of the incident radiation 
for the working standards and sphere coatings used. It follows 



that the random error in F c j is given by: 



8F r 



.p 



\Qc 



+ 



fifty 



+ 



8F, 

F 



?)T- 



(42) 



The largest new contribution to the uncertainty which is 
introduced in this step is 8Qj. This is due to the variability 
between the representative samples and the resulting uncer- 
tainty in whether these samples properly represent the sphere 
coating. If the values of Q c and Qj differ significantly a small 
correction for instrument non-linearity must be made in the 
ratio QclQj- With this non-linearity is associated a small 
contribution to the systematic uncertainty. 



AF = F • 

LX1 C,J l C,3 



IA (Q c /Q } ) ) 



+ (9j 



(43) 



Three different working standards are currently used. 
These are identified and described in appendix D. The values 
of F c j for each of these standards for each of three determina- 
tions are given in tables XI, XII and XIII, and sample 
calculations for the uncertainties appear in table XIV. 
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Table XI 



8 

A 

€ 
W 

W 



2.2 x 10- 

3.3 x 10- 
2.2 x 10" 
4.5 



Reflectance Factor of Working Standard 1 
(Russian Opal Glass, Polished) 



Wavelength 
(nm) 


F U1 


Fl.2 


^1,3 


Fi* 


*u 


400 


0.9793 


0.9734 


0.9738 


0.9759 


0.9753 


425 


0.9781 


0.9739 


0.9743 


0.9756 


0.9753 


450 


0.9803 


0.9765 


0.9779 


0.9781 


0.9781 


475 


0.9841 


0.9807 


0.9812 


0.9822 


0.9819 


500 


0.9854 


0.9823 


0.9846 


0.9836 


0.9838 


525 


0.9856 


0.9821 


0.9845 


0.9836 


0.9838 


550 


0.9837 


0.9813 


0.9830 


0.9823 


0.9835 


575 


0.9810 


0.9785 


0.9805 


0.9796 


0.9798 


600 


0.9792 


0.9761 


0.9785 


0.9775 


0.9777 


625 


0.9777 


0.9742 


0.9769 


0.9757 


0.9760 


650 


0.9770 


0.9745 


0.9774 


0.9756 


0.9761 


675 


0.9771 


0.9751 


0.9772 


0.9759 


0.9763 


700 


0.9761 


0.9736 


0.9761 


0.9747 


0.9750 


725 


0.9742 


0.9718 


0.9741 


0.9728 


0.9732 


750 


0.9716 


0.9687 


0.9711 


0.9699 


0.9702 



1.6 x 10- 
5.5 x 10" 

1.7 x 10" 
6.0 



2.7 X 10" 

3.3 x 10" 

2.7 x 10" 
3.7 



1.3 x 10- 
5.5 x 10~ 

1.4 x 10" 

10.5 



1.2 x 10- 
5.5 x 10" 

1.3 X 10" 

14.2 



Table XII 







Reflectance 


Factor of Working Standard 2 












(Vitrolite) 








Wavelength 
(nm) 


F2.1 


^2,2 




^2,3 


^2,2 


F 2 . 3 


400 


0.9100 


0.9049 


0.9049 


0.9071 


0.9065 


425 


.8944 


.8901 




.8907 


.8920 


.8916 


450 


.9006 


.8977 




.8981 


.8990 


.8987 


475 


.9140 


.9116 




.9083 


.9126 


.9114 


500 


.9192 


.9169 




.9178 


.9179 


.9179 


525 


.9242 


.9204 




.9226 


.9221 


.9222 


550 


.9245 


.9222 




.9243 


.9232 


.9235 


575 


.9215 


.9190 




.9213 


.9201 


.9204 


600 


.9160 


.9125 




.9155 


.9140 


.9145 


625 


.9099 


.9066 




.9085 


.9080 


.9082 


650 


.9050 


.9027 




.9043 


.9037 


.9039 


675 


.9035 


.9013 




.9040 


.9023 


.9028 


700 


.9010 


.8991 




.8999 


.8999 


.8999 


725 


.8959 


.8939 




.8953 


.8948 


.8949 


750 


.8893 


.8867 




.8882 


.8878 


.8879 


8 


2.2 X 10- 3 


1.6 X 10- 3 


2.3 


X 10- 3 


1.3 X 10- 3 


9.8 X 10- 4 


A 


3.2 X 10- 4 


5.3 X 10- 4 


3.2 


x 10- 4 


5.3 X 10- 4 


5.3 x 10- 4 


e 


2.2 X 10- 3 


1.7 X IO- 3 


2.4 


x 10- 3 


1.4 x 10" 3 


1.1 x 10- 3 


w 


4.5 


5.9 


4.2 




— 


— 



w 



10.4 



14.6 



The systematic uncertainty associated with the F c for a 
particular standard is, by definition, independent of the 
number of determinations which are made. However, the 
random uncertainty associated with a given F c can be re- 
duced by taking the average of the F c 's from a number_ of 
determinations. We maintain a running weighted average F c j 
of the results of all determinations from 1 through y by means 
of the calculation 



where the weighting factor w c 
F c is 



for the /th determination of 



= 1/100 eF c 



(45) 



where eF Ct j is the total uncertainty in F Cy j as is given in 
appendix B. The weighting factor W c j for the average of they 
determinations is 



Fcj = 



^C,j-M Cj— 1 ' W C,j^ C,j 

~w7~ 



(44) 



W cj ~ X w ci ■ 



(46) 



45 



Table XIII 

Reflectance Factor of Working Standard 3 
(Porcelain Enamel 80-1) 



Wavelength 
(ran) 



3,2 



3,3 



400 


0.7615 


425 


.7936 


450 


.8042 


475 


.8093 


500 


.8105 


525 


.8104 


550 


.8090 


575 


.8063 


600 


.8030 


625 


.7990 


650 


.7949 


675 


.7926 


700 


.7909 


725 


.7885 


750 


.7855 


8 


2.1 X 10 


A 


3.0 X 10" 


e 


2.2 X 10 


w 


4.6 


W 


- 



0.7573 


0.7571 


.7911 


.7915 


.8019 


.8034 


.8072 


.8073 


.8078 


.8101 


.8091 


.8096 


.8083 


.8094 


.8045 


.8071 


.8017 


.8032 


.7969 


.7994 


.7937 


.7960 


.7911 


.7925 


.7900 


.7927 


.7869 


.7886 


.7839 


.7860 


1.5 x 10~ 3 


2.3 X 10" 


4.7 X 10~ 4 


3.0 X HT 


1.6 X 10- 3 


2.3 X 10" 


6.2 


4.2 



0.7591 
.7922 
.8029 
.8081 
.8090 
.8096 
.8086 
.8053 
.8022 
.7978 
.7942 
.7917 
.7904 
.7875 
.7846 



0.7585 
.7920 
.8030 
.8079 
.8093 
.8096 



.8058 
.8025 
.7983 
.7947 
.7919 
.7910 
.7878 
.7850 



1.3 x 10- 
4.7 X 10" 

1.4 X 10" 

10.8 



1.1 x 10" 
4.7 X 10- 

1.2 x 10- 

15.0 



Table XIV 
Calculation of the Uncertainties in the p(6°,A.) to F(6°,k) and F(6°,A.) to F c , j(6°, A.) steps as determined for Working Standard No. 2. 





Sym 


bol 




Value 




i 




BaS0 4 [14] 


Halon [15] 


BaS0 4 [16] 




(a) 


8F(6°, X) 
(-op(6 , a)) 


9.5 X 10- 4 


7.9 X 10- 4 


2.1 X lO" 3 


See table X, entry (f). 


(b) 


Fi 




0.9767 


0.9932 


0.9806 


See tables IV, V, and VI. 


(c) 


dQ c 




2.7 X 10~ 4 


8.7 X lO" 4 


3 X lO -4 


From three measurements in each determina- 
tion. 


(d) 


Qc 




0.929833 


0.931036 


0.92922 


Measured. 


(e) 


SQi 




2 X 10- 3 


1.2 X lO" 3 


1.2 x lO" 3 


From measurements on a number of samples. 


(0 
(g) 


F2,J 




0.985022 
.9245 


1.00283 
0.9222 


0.985982 
0.9243 


Measured. 

Calculated from (b), (d) and (f) above in equa- 
tion (41). 


(h) 


SFt.i 




2.2 x 10-3 


1.6 x lO" 3 


2.3 X lO" 3 


Calculated from equation (42) and (a) through 
(g) above. 


(i) 


HQe/Qi). 

(~- A<?i) 


-1 x 10- 4 


-1 x lO" 4 


-1 x lO -4 


From reference [17]. 


(j) 


AF, 




3.3 x 10" 4 


5.6 X 10" 4 


3.3 x 10" 4 


From tables IV, V, and VI. 


(k) 


A^2, i 




3.2 x 10- 4 


5.3 X lO -4 


3.2 X 10" 4 


Calculate from equation (43) and (b), (d), (f), 
(i) and (j) above. 


(1) 
(m) 
(n) 
(o) 

(P) 

(q) 
(r) 


^2, j 
W c , j 
Vc.3-1 

F«s 




2.2 X 10-3 
4.5 


1.7 x lO" 3 
5.9 
10.4 
0.9232 

1.3 x lO" 3 

5.3 X lO -4 

1.4 x 10" 3 


2.4 X 10" 3 
4.2 
14.6 
0.9235 

9.8 X lO -4 

5.3 X 10" 4 
1.1 x lO" 3 


Quadrature combination of (i) and (k) above. 
Calculate from equation (45) using (1) above. 
Calculate from equation (46) using (m) above. 
Calculate from equation (44) using (g) and (m) 

above. 
Calculate from equation (47) using (g), (m), 

(n), and (o) above. 
Largest entry from (k) above. 
Quadrature combination of (p) and (q) above. 
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The random uncertainty in the F c j is given by 



8F r , = 



[(W^SFcj-,) 2 + (w CJ 8F cj yr* 



Wr 



(47) 



The values ofF Cf2 andF Ct s appear in tables XI, XII, and XIII 
along with the F c j. The results of these three determinations 
are also displayed by way of summary in figure 16. 

V. Conclusions 

A. Precautions and Pertinent Parameters 

In the work reported in this paper, a large number of 
potential sources of error in determining reflectance factor 
using the Van den Akker method were investigated in order 
to evaluate the magnitude of their importance. Many of these 
sources of uncertainty were of little importance in themselves 
and even cumulatively were of importance only to a labora- 
tory seeking the highest accuracy possible. For a laboratory 
for which a total uncertainty on the order of ±0.005 is 
tolerable, a number of the detailed steps described in section 
IV are not necessary. In this section we will briefly describe 
precautions which should still be taken in the case of such 
truncated measurement procedures and make suggestions for 
possible improvements in the overall procedure. 

1. Determining the Van den Akker Reflectance p v (A.) 

In an experiment designed and carried out with reasonable 
care, the Van den Akker reflectance p v can be determined 
very accurately. The principal strong point in this method is 
the insensitivity of the basic measurement accuracy to uncer- 
tainties in any of the measured parameters such as the port 
diameter, the sphere diameter, or the ratio of the reflectance 
of the sphere to the reflectance of the wall sample. The 
sensitivity to such parameters is kept lowest by using a large 
sphere with a small port (low f) and by using a highly 
reflecting coating in the sphere so that the sphere reflectance 
is still reasonably high. 

Considerable care should be taken however to make the 
sphere coating in such a way that the Van den Akker reflec- 
tance can be related meaningfully to the reflectance of a 
sample of the coating. This means that a coating should be 
chosen the reflectance of which is uniform and very repro- 
ducible. It is especially important that no cracks appear in 
the finished coatings and that the coating should be thick 
enough to be opaque. In this regard, it is probably best not to 
incorporate a removable port in the Van den Akker sphere at 
all, but rather to rely upon the reproducibility of the coating 
and use the average Q from several separate samples of 
coating for Q t in eq (6). 

2. Adjustment from the Van den Akker Reflectance p r to the d/h 
Reflectance p(6°, X) 

As can be seen from figure 13, the magnitude of this 
adjustment was in no case greater than —0.006. Since p(T, 
X) can never be greater than 1, it follows that this adjustment 
will be smaller for samples for which p v closer to 1. It also is 
clear that the need for this adjustment comes about from the 
rise in reflectance at near-grazing incidence. This rise is due 
to the specular reflectance of the slightly glossy surface which 
is produced in pressing. For this reason, a coating technique 



which produces a rougher or more matted coating is to be 
preferred if no adjustment is to be made. It is probably the 
roughness of the coating that Goebel, et al., [2] produced by 
scraping which led to the rather remarkably good agreement 
between their results, which are the basis of the earlier NBS 
scale of F(6°, A.), and the present adjusted results. (See fig. 
17 and the discussion in section V.C. which follows.) An 
error of as much as 0.01 in the upward direction can result 
from assuming p v and p to be equal, and an adjustment based 
on relative d/h reflectance measurements NJ(l, A.) should be 
made if an uncertainty of less than ±.01 is desired. In all but 
the highest accuracy work, it is safe to assume the self- 
radiance is Lambertian, i.e., V^y, X)/cos y is constant. 

It should be pointed out that in every case p v and p^n 
are very nearly equal, as can be determined from the C(T, k) 
curves in figure 12. Therefore the Van den Akker method is 
especially well adapted for calibrating measurements of di- 
rectional-hemispherical reflectance at 45° incidence. 

3. Adjustment from d/h Reflectance p(6°, A) to d/h Reflectance Factor 
F(6°,X) 

Since most reflectometers compare reflectance factor F(6°, 
X) rather than reflectance p(6°, X), it is in principle necessary 
to make an adjustment. However, since the solid angle 
subtended at the sample is small compared to the total 
hemisphere, any departure from Lambertian reflectance by 
the sphere coating would have to be large if the difference 
between p and F is to be significant. In the sphere coating 
materials ordinarily in use, the reflectance factor departs 
greatly from 1 only for small solid angles of collection at very 
small observation angles /3 (fig. 14). Therefore in general p is 
less than 0.001 higher than F and this adjustment can be 
ignored or estimated from data in the literature [13]. 

Note that if the object of the measurement is to determine 
the radiative transfer properties of the sphere coating itself, 
then p is the quantity of interest and the adjustment to F 
should not be made. This would be the case, for example, if 
the coating were being studied in order to determine the 
radiative heat transfer to the coating material. 

4. Determining the d/h Reflectance Factor F(6 ", A) of Working 
Standards 

Next to the adjustment from p v to p, it is this step which 
contains the highest potential for error. These measurements 
are meaningful only to the extent that the samples of coating 
with which the working standards are compared are repre- 
sentative of the sphere wall coating. It is for this reason that 
the sphere coating and the samples should be prepared in 
precisely the same way and that reproducibility of reflectance 
is an important requirement of the coating. It is also this 
reason that favors Van den Akker's original technique of 
having many removable sections in the sphere wall itself, 
provided that the presence of such removable sections does 
not give rise to irregularities in the sphere coating. 

B. The Present NBS Scale of d/h Reflectance 
Factor F(6°, X) 

The present NBS scale of spectral d/h reflectance factor is 
maintained by means of three carefully preserved working 
standards which have been evaluated at 25 nm intervals. The 
reflectance factor data for the first three determinations are 
given in tables XI, XII and XIII. Figure 16 provides a 
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WORKING STANDARD NO. I 
(RUSSIAN OPAL GLASS) 




• l»t RUN (BaS0 4 ) 

X 2 nd RUN (HALON) 

O 3rd run (BaS0 4 ) 

S*^^ AVERAGE |*t&2 ndRUN 

• WORKING STANDARD NO. 2 
(VITROLITE) 




400 



500 600 
WAVELENGTH (nm) 



700 750 



FIGURE 16. Reflectance factor differences for the measurement of three 
working standards in three reflectance scale determinations. 



graphical record of the establishment of the average value for 
the three determinations. The representative error bars given 
to the left of the diagram actually were derived for 550 nm, 
but are fairly representative of the entire spectral range. The 
zero line in each case represents the average for the three 
determinations, and the shaded area represents the range of 
the total uncertainty of that average. Although the distribu- 
tion of the points in general is consistent with the uncertainty 
levels which have been assigned, there is a lack of consist- 
ency from standard to standard between the relative values 
from pairs of determinations. For example, for working stan- 
dard No. 1 (Russian opal glass), the third determination 
reflectance values are consistently lower than those of the 
first determination, whereas, for the other two working stan- 
dards the first and third determinations yielded values which 
were more nearly equal (fig. 16). This difference in general 
behavior is too large to be explained by the uncertainties in 
the measured values of Q c in expression (41). The most 
straightforward explanation for this difference is a slight drop 
in the reflectance of the working standard No. 1 during the 
interval between the second and third determinations. A 
better knowledge of the stability of the working standards will 
be obtained as more data are gathered from additional deter- 
minations. 

The results of the work reported in this paper document the 



establishment of a scale of d/h reflectance factor at 25 nm 
intervals over the wavelength range 400 to 750 nm. The error 
analysis indicates that for high quality, uniform samples with 
reflectance greater than 0.5 there is only a very small proba- 
bility that our measured reflectance factor values will be in 
error by more than 0.0015. 

C. Relationship to Other Scales 

In order to determine the relationship between the newly 
established NBS scale of d/h reflectance factor and the 
former one, we measured the three working standards on the 
NBS scale established in 1965. The results of this compari- 
son are summarized in figure 17. In this figure, the data 
points indicate the departure of the old scale from the new for 
each working standard. The departure of the average value of 
the spectral 6°//i reflectance factors for the three standards on 
the 1965 scale from the same average on the new scale is 
shown by the light continuous line (GE average). The shaded 
area in the figure represents the uncertainty of the new scale, 
and the error bar on a central point indicates the uncertainty 
attributed to the 1965 scale. On the average, the old scale 
appears to depart from the new by about 0.002 at the short 
wavelength end of the range with the magnitude of the 
departure diminishing as the wavelength increases. 



.005 




uj -.005 
o 

z 

UJ 



-.01 



-015 






x x X x 

X 



X X 
""X" 



• RUSSIAN OPAL GLASS 
X VITROLITE 
H- ENAMEL 80-1 
"\^ GE AVERAGE 




NRC VALUE 



400 



500 600 

WAVELENGTH ( nm ) 



700 750 



FIGURE 17. Comparison of present NBS scale of F 6 o_ h (taken as zero line) 
to previous NBS scale and to Results of an intercomparison with NRC. 



In 1974, measurements were made of several samples of 
pressed BaS0 4 and matte ground Russian opal glass at the 
National Research Council Laboratories (NRC) in Canada 
and at NBS on the 1965 NBS scale. The heavy line in figure 
17 (NRC value) departs from the light line (GE average) by 
the difference between the averages of the reflectance values 
determined in the intercomparison. In this way, the present 
NBS scale and the NRC scale of 1974 are compared indi- 
rectly. The rather noisy appearance of the difference can very 
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probably be attributed to the impreeision of the NBS GE 
spectrophotometer. 

The clashed line at approximately —0.004 indicates the 
combined error for this comparison. Since the NRC value 
differs from the present NBS value by over twice that amount, 
it is clear that at least one of these scales is in error by an 
amount greater than that claimed for it. At the time of this 
writing, work is under way on a new direct intercomparison 
between the present NRC and NBS scales of directional- 
hemispherical reflectance factor and on joint efforts to deter- 
mine the cause of any differences. 

D. Looking Ahead 

We have already begun work on a determination of d/h 
reflectance factor by the Sharp-Little method [25, 26]. Most 
of the data have been obtained, and detailed analysis will be 
begun soon. Preliminary results indicate that the reflectance 
factors of the three working standards obtained by that 
method agree with those obtained by the Van den Akker 
method to well within the uncertainty of 0.0015. Therefore, it 
appears that the present NBS scale is accurate to within the 
uncertainty reported in this paper. 

Mainly because of the slight instability and non-uniformity 
of the working standards, it is doubtful that an uncertainty in 
F(6°, A.) of less than 0. 1 percent of its value can be achieved 
on a practical, routine basis by any technique which relies on 
such standards. However, these working standards are more 
stable and more uniform than most samples encountered in 
practice. For this reason, we believe the newly established 
NBS scale of d/h reflectance factor to be sufficiently accurate 
and precise to serve the current needs for such measurements 
in science and industry. 
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VII. Appendix A. Accounting for 
Retroreflectance 

Retroreflectance refers to the tendency of objects to prefer- 
entially reflect radiation hack in the direction from which it is 
incident. This can arise from several causes. If a focussing 
element such as a small dielectric sphere directs radiation 
onto a reflecting surface at or near its focal point, the 
reflected radiation will he returned back into the direction 
from which it came. This type of retroreflectance can be very 
directional, i.e., most of the radiation returns within a few 
tenths of a degree of the angle from which it came. A second, 
much less directional retroreflectance effect, results from the 
fact that shadows are not visible from the direction of inci- 
dence. This mechanism is important only for surfaces which 
are not highly reflecting, so that the radiation is not re- 
reflected out of the shadows. In the case of the highly 
reflecting coatings being used in the spheres in this experi- 
ment, the highly directional retroreflectance dominates. 

7 . The Effect of Retroreflectance 

In order to obtain an estimate of how the presence of 
retroreflectance will affect a determination of the Van den 
Akker reflectance p v , we will use a simple model in which 
most of the radiation is reflected in a Lambertian distribution 
but a small retroreflected fraction f r appears above the Lam- 
bertain background. The effect which the retroreflectance has 
upon the flux </>i(X) (see III. A.) striking the sphere boundary 
can be seen by going one step backwards in determining the 
source of fa. The flux fa comes from the walls by reflection, 
and the incident radiation for this reflection is from the 
sphere walls, but not from the port. Therefore, in the pres- 
ence of retroreflection the irradiance on the sphere boundary 
is higher on the wall coating than in the port. The makeup of 
fa can be determined by noting that, in this model, the flux 
incident upon the walls (which also comes from the walls) is 
either absorbed or reflected. The total reflected flux, which is 
fa , is made up of two parts 



fa - fa + </>i 



(Al) 
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where (f) r is the retroreflected fraction of the reflected flux 

4> r =Mn (A2) 

and fa is the diffusely reflected fraction 



<t>'i = (1 -fr)<l>i' 



(A3) 



With this model, the power balance equation, equation (3) of 
part III. A., takes on the following form. 

Pv( 1 +/rp;)(l - f)pt4>0 (A4) 

= ^+[(l-M + f](l- p' v ) 

where p t is reflectance of the target exclusive of the aug- 
mented reflection in the retro direction, which is lost through 
the port of the measuring instrument in every case. The 
relative reflectance Q s of the sphere is 

Q, = kWj{\ -f) + 4> (l + frP'MH ~ /')] (A5) 

and the relative reflectance of the target is, as before 

Qt = k<f> oPt (l-f'). (A6) 

Equations (A2), (A3), (A4), (A5), and (A6) can be solved to 
obtain an expression for p v in the following form: 



Pv 



1 -/ (1 -f't _l_ 

1 "/(I ~fr) ' 



a 



/( l + P'yfrtf "/"I 
(1-/) 



1 - a 



#r(l -frPv 2 ) 



(A7) 



1 - P i[l -Al -fr)])Qs 

(A7a) 



Comparing (A7) with (6), it is seen that fin the dominant first 
term is replaced by/(l — f r ). This in essence states that the 
fraction the flux escaping from the port has been reduced 
relative to the flux lost to the walls by an amount proportional 
to the retroreflected fraction. The correction term a for the 
singular treatment of the first target reflection when the target 
is in the auxiliary sphere is also modified by the retroreflec- 
tance, making the form of this term so complicated that a 
closed solution in terms of the measurement parameters is not 
practical. We solve the equation iteratively, using as a first 
approximation for p' v the value obtained from (A7) with a = 
0. Taking the retroreflectance into account has a very small 
effect on the calculated value of p v . For example, if/ r = 
0.001, the effect of not ignoring this in the case of a sample 
for which p' v —0.98 is to reduce the calculated p' v by 
approximately 10 -6 . Note that p' v must be the total reflec- 
tance, including retroreflectance, since 1 — p' v represents 
loss of radiation by absorption. 

2. Models for Retroreflective Mechanisms 

Our measurements of retroreflective bidirectional reflec- 
tance factor, the data for which are shown in figure 15, were 



limited to observation angles greater than 0. 1°. Therefore we 
sought a model which could be used to interpolate our data to 
0° observation angle so that we could evaluate the integral in 
eq (39). 

A model based on shadowing was developed by Hapke [27] 
to explain the retroreflective phenomena observed from the 
lunar surface. This model, which was quite successful in 
Hapke' s application, was used by Egan and Hilgeman [13] in 
an effort to quantitatively evaluate the retroreflectance which 
they observed from barium sulfate paint and other highly 
reflective white coatings. Other possible mechanisms for 
producing retroreflectance include cube corner reflectors 
formed from broken cubic crystals and retroreflectance by 
focussing, as occurs in glass beaded retroreflecting sheeting. 
In the following paragraphs we will treat special cases of each 
mechanism. 

a. Shadowing 

The Hapke model includes as one of its basic assumptions 
that the absorbance of the scattering particles is very high. In 
this way, there will be a considerable contrast between the 
radiance coming from shadowed areas in which the radiation 
undergoes several reflections and the radiance reflected from 
an unshadowed area. 'Although there is relatively little loss in 
highly reflecting materials such as barium sulfate or Halon, it 
is still possible for an initial first surface reflection to add to 
the background of multiply reflected radiation in a preferen- 
tial manner. In order to estimate the order of magnitude of the 
non-diffuse reflection, one can use layers of close-packed 
spheres as a model for the many randomly oriented particles 
in the coating. Referring to one such sphere as shown in 
figure 18a, one can calculate the reflectance factor for such 
an array as follows. The bidirectional reflectance factor is 
given by 



F(V, u) = 



ttL(u) 



(A8) 



where E(U) is collimated irradiance incident in the direction 
U and L(u) is the reflected radiance in the direction u, both 
averaged over a sufficiently large area of the sample. In the 
case of retroreflectance, we set the unit vectors U and u 
equal. The average reflected radiance is given by definition 



i(u) = 



3>(o>, A) 
o)(u)Acosy ' 



(A9) 



where co(u) is an element of solid angle oriented in the 
direction (u), A is the area of the surface over which the 
average is being taken, and *&((*), A) is the flux reflected from 
area 5 into solid angle (t). Collimated incident radiation will 
come to a virtual focus at a point / which is r/2 from the outer 
surface of the sphere. The solid angle 0) is defined in terms of 
a small arbitrary area a on the sphere 



(O 



4a 

,-2 



(A10) 



The flux reflected from this area into (O is given by 

$ = pEa/ cosy (All) 
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where p is the reflectance of the sphere surface. (Note that 
since we are treating retroreflectance, a is chosen so that it is 
centered on a radius in the direction of E, and therefore the 
direction of incidence is normal to a.) The area A occupied by 
this sphere in a hexagonal close-packed array is 



A = 



6r 2 



VS" 



(A12) 



Using expressions (A9), (A10), (All), and (A12) above in 
(A8) one obtains under conditions of retroreflectance: 

F (U, u) = — p/cos^y = 0.227 p\ 'cos 2 y. (A 13) 
24 

This would be roughly the maximum amount of additional 
bidirectional reflectance factor which might occur in the 
retro-direction under ideal conditions of shading by sur- 
rounding particles. The cos 2 ^ term in the denominator is an 
artifact of the "array of spheres" model we are using. In the 
continuous surface being represented by the array of spheres 
model, the projected surface area goes as cos 2 y so that the 
cos 2 *)/ term should be omitted. 

b. Cube corners 

In the ideal case, the bidirectional reflectance factor for a 
cube corner reflector would be infinite at zero observation 
angle and zero at all other observation angles. However, 
imperfections in the cube corners and diffraction tend to 
spread the reflected radiation out so that bidirectional reflec- 
tance factor has meaning in describing reflection by cube 
corners. For the model for maximum reflectance, we will 
assume that the cubes are perfect and that diffraction is the 
limiting factor. For near normal incidence the total retrore- 
flectance from a close-packed array of cube corner reflectors 
is 



p t = p?(l - p 2 f 



(A14) 



where pi is the reflectance of the cube faces and p 2 is the 
reflectance of the entrance surface. For a perfect internal 
corner pj is 1 and p 2 is roughly 0.05, so that p t is roughly 
0.9. For an external corner (cube corner "hole"), p 2 — and 
p x is roughly 0.05 so that p t is roughly 1.25 X 10 -4 in that 
case. The reflected radiance is distributed in the diffraction 
pattern associated with the reflecting faces treated as aper- 
tures. For this order of magnitude calculation, we will assume 
the pattern to be that of a circular aperture of effective radius 

[A 
r = \i — where A is the area of the entrance surface. The 

* IT 

radiance in such a pattern would be given by [28]. 



/ x EA 



2Jt(x) 



(A15) 



x = 2tt r sin (;8)/A, 
from which il follows that the reflectance factor is given by 



F(j3) = 



irA 



2/,(*)> 



x J 



(A16) 



The data in figure 15 suggest that the particle sizes are such 
that the central bright spot would subtend an angle of about 1 
for BaSCV The edge of the central bright spot falls at x = 
3.8. For an effective wavelength of approximately 550 nm, it 
follows that the particle size would have to be such that r — 
20 ptm. Thus, the reflectance factor at /3 < 1° would be 
roughly 



F(/3) - 10V 



(A17) 



Focussing retroreflectors 



If the index of refraction of a bead is such that the radiation 
passing through it comes to a focus behind it, a suitable 
reflector may be placed at the focal point (I in fig. 18c.) and 
the radiation will be reflected back through the bead into a 
collimated retroreflected beam. Here, the spreading of the 
reflected radiance can be brought about either by diffraction, 
poor focussing, or both. The diffraction limiting case in the 
previous paragraph can be applied to the Halon data in figure 
15. In this case, the central bright spot appears to subtend an 
angle of about 0.5° which corresponds to an effective lens 
radius of r —40 p,m. If one assumes a refractive index of 1.5 
so that the radiation will come to a focus on the back of the 
bead and from experience with cheap cameras assumes that 
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C. 
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FIGURE 18. Diagrams for estimating magnitude of retroreflectance by 
various means. 
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the diameter of the effective lens area is roughly one-eighth of 
the diameter of the bead, it follows from expressions (A 12) 
and (A 16) that a closely packed array of such spheres would 
have a reflectance factor for /3 <^ 1 of 

F(j3) - 6 X 10 2 Pt (A18) 

for a specular reflector of reflectance Pt located at I and 

F(J3) - lOp, (A19) 

for a diffuse reflector of reflectance Pt at I. 

3. Interpolating the retroreflectance data 

From the foregoing analysis, one can make the following 
deductions, assuming that the 0.5 excess reflectance factor is 
caused by only one mechanism: 

a. Shadowing cannot by itself account for the observed 
retroreflectance . 

b. At least 6 X 10 -5 of the area is effectively taken up 
with inside corner cube retro reflectors. 

c. At least 0.4 of the area is effectively taken up with 
outside corner cube retro reflectors. 

d. At least 1.6 X 10 -2 of the area is effectively taken up 
with focusing sphere retroreflectors backed by specular 
reflectors of reflectance 0.05. 

e. At least 5 X 10~ 2 of the area is effectively taken up 
with sphere retroreflectors backed by diffuse reflectors 
of reflectance 1. 

Since the mechanisms are not mutually exclusive, it is ex- 
pected that all will operate to some extent in pressed powder 
samples. However it is highly unlikely that cube corner 
reflectors will be formed in Halon, so it would be expected 
that the mechanism in d. or e. above would predominate. 
Even though the cleavage of a crystal in such a way that it 
forms an inside corner cube reflector is not very likely, 
mechanism b. probably predominates in the BaS0 4 powder. 
Shadowing will contribute some effective retroreflectance in 
both cases. 

The predominate mechanisms appear to depend on the 
sum of the diffraction patterns associated with a number of 
reflectors, each with a different effective aperture, and the 
statistical distribution of the aperture sizes is unknown. 
Therefore we chose to use an emperical fitting function. The 
function 

F(/3) = 1 + exp(a 4- b/3 c ) (A20) 

was found to fit the data adequately. Fitting all BaS0 4 data, 
we obtained 

a = -0.393 
b = -0.966 
c = 0.57 

and fitting the first seven points of the Halon data we ob- 
tained 

a = 1.18 
b = -3.41 
c = 0.21 



These fitting functions are shown as solid curves in figure 15. 
The average of the function F(/3) for BaS0 4 over the range 
< /3 ^ 0.5° is 1.42 which is in good agreement with the 
data of Egan and Hilgeman [13]. The difference values for 
p(6°, k) — ^(6°, k) were obtained by using the fitted func- 
tions in eqs (39) and (40) with r = 1.9 cm and a 1 ' = 29 A cm. 
In each case, the contribution from the region of extrapolation 
(|j3| < 0.1°) was only a small part of the total flux in the 
retroreflected component, so the choice of fitting functions 
was not critical. A crude measurement of BaS0 4 using a 
beam splitter yielded p(6°, k) - F(6°, k) = 0.0006 which is 
in good agreement with the value of 0.0007 obtained from the 
data in figure 15. 



VIII. Appendix B. Interpretation of Error 
Statements 

The error analysis in this paper involves three distinct 
classes of errors. The first class of error is the measured 
random error. This type of error appears as noise in the 
measurements and is evaluated by making each measurement 
several times and performing the usual error analysis on the 
results. For this type of error, we will cite three times the 
standard deviation of the mean. 

The second class of error is the estimated random error. 
This class of error contributes to the outcome in the same way 
as the measured random error, but, either because of the 
inconsequential effects of the error or because of experimen- 
tal difficulties which would be involved in evaluating it, we 
choose to estimate the size of the error instead of determining 
it by making several measurements. In making such an 
estimate, we choose to estimate the smallest range of uncer- 
tainty within which the correct value will fall with a very high 
probability, i.e., an intuitive ninety-five times out of a 
hundred. For this type of error, we will cite the size of the 
estimated error itself. The estimation process is carried out in 
such a way that the sign of the error is indeterminant. If it is 
known that a particular cause of error will result in an error 
which is always of one sign, we will correct the measured 
value to a most probable value and reduce the magnitude of 
the estimated random error appropriately. The estimated 
random error is combined in quadrature with the measured 
random error to obtain the total random error. All random 
errors, whether measured or estimated, will be denoted by 
the symbol 8. 

The third type of error is the estimated systematic error, 
indicated by the symbol A. This error is like the estimated 
random error in all respectes except one, namely it is not 
independent from determination to determination, i.e., it will 
affect all determinations in the same way. It should be 
pointed out that this does not mean that the sign of the 
measurement is known, but only that the sign of the effect is 
known to be fixed from determination to determination. For 
example, in correcting for the reflection of radiation from the 
beveled lip of the sphere, we estimate that the radiation will 
be reflected with an average reflectance intermediate be- 
tween one and the normal reflectance of the stainless steel. 
(See discussion preceding (19).) If the actual reflectance is 
less than our estimated value, determinations of reflectance 
made using this estimate will result in p v values which are all 
slightly higher than if the estimates were correctly made. In 
the course of many determinations, the total random error will 
be reduced, but the systematic error will remain. Since the 
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sign of the systematic errors is not known, systematic errors 
from independent causes will be added in quadrature, and at 
each point in the analysis at which it is desirable to estimate a 
total uncertainty, the systematic error and the total random 
error will be added in quadrature. However, when proceeding 
to combine the results of several determinations, we will first 
combine the random errors in the ordinary way and add the 
largest systematic error in quadrature at the end. In this way, 
the calculated systematic error is not reduced by repeating 
the measurements. The total uncertainties so calculated are 
to be regarded as determining a range within which we expect 
the correct value to fall with a probability of approximately 
0.95. This total uncertainty will be indicated by the symbol 



IX. Appendix C. Calculating C(r, k) 

The integrals in (11) are evaluated in closed form after 
fitting the experimental data for N g and N G . 



In the case of the original knife-edge data, the values of the 
coefficients, as obtained by least squares fitting, are given in 
table CII as describing the "upper limit" curve. In the 
original data, the least uncertainty was associated with the 
r = reading. Therefore, all data were normalized to the T 
= reading and b was taken to be 1 when the data were 
fitted. Because of the symmetry in azimuth, b\ and b s are 
taken to be zero. 

In the accepted data, N g (T)/ cos T was taken as 1 up to F — 
7r/4. To obtain the remaining part of this curve we obtained 
N g (y) according to the expression 



w;(y)/cosy= 1 + dy. 



(C3) 



The difference data N g (y) ~ N g (y) using the data for the 
angles 7r/3 (60°), 5tt/12(75°) and 177t/36(85°) were fitted 
using the form 



n\, 



ty,(y)/(cos y) = 1 + b'iT - -r + b' 2 (T 

4 



,) 4 . (C4) 



7. Fitting N G (T, X) 

The data for N G were fitted by an expression of the form: 

7V G (r, X) = a (X) + ChMP + a 2 (k)F\ (CI) 3. Calculating C(I\ k) 



After b 1 and b 2 were determined, these were used to obtain 
the bi in the expression (C2). 



An even function of T was used, since we have assumed the 
directional-hemispherical reflectance of the coatings does not 
depend on the angle of azimuth of the incident direction, but 
depends only on the angle of elevation. The coefficients in 
(CI), as determined by a least squares fitting of the experi- 
mental data are given in table CI. 

2. Fitting N (T, X) 

Since the experimental data for N g (y, X)/cos (y) were found 
to be nearly independent of wavelength, only one function 
N g (T) is needed for each type of coating. However, since 
there was such a large uncertainty in the N g (T) data, three 
sets of N g (T) were fitted as described in part IV. B. 3. in the 
main text. The form used for the fitting the data was: 



With the expressions (CI) for /V G (I\ A) and (C2) for N g (T, 
X) substituted into eq (11), the integrals in that equation can 
be evaluated in closed form in terms of integrals of the form 



h = Jo /2 r* cos r sin r<*r 

for the upper and lower limit functions, and 
/« = Jo' 4 P cos r sin TdT 
and 



ir=j?£r<cosrsinr<fr 



(C5) 



(C6) 



(C7) 



a^d/cost 



1=0 



(C2) 



Table CI 

Coefficients for the expansion of /V^T, A.) in even powers of the angle of 
incidence T in radians at three different wavelengths for three sphere wall 
coatings. 



Coating 



BaS0 4 [14] 



Halon [15] 



BaS() 4 [16] 



Wave- 
length 


a 


X 10 3 


«2* 

x 10 3 


450 
550 
750 


0.999942 

.99978 

1.00011 


5.99 
4.40 
3.27 


0.438 

0.159 

-0.00814 


450 
550 
750 


0.999875 

1.00004 

1.00027 


3.97 
4.17 
3.68 


-0.172 
-0.985 
-1.20 


450 
550 
750 


0.999829 

1.00047 

1.00012 


4.08 
1.43 
3.91 


0.405 

0.741 

-0.471 



* The second and third digits are not significant but are kept to avoid 
round-off error in further calculations. 



for the accepted value functions. The values of theseintegrals 
are given in table CIII, and the resulting values of JS G (\) are 
given in table CIV. It can be seen from table CIV that the 
difference between the various assumptions concerning the 
self- radiance has a very small effect on the value of N G in all 
cases. Therefore only the accepted value oiN G will be used in 
each case to calculate C(T , A.). The results of the scanning 
which has been done to date with the telescopic detector 
indicate that the real value probably lies toward the lower 
limit from the accepted value. On this belief, we will assign 
the largest difference between an accepted value and a lower 
limit value, 0.00005, as the magnitude of the uncertainty in 
Nq due to the uncertainty in evaluating the self-radiance by 
measuring N g . 

X. Appendix D. The Three Working Standards 

The properties which standards for reflectance factor mea- 
surements should have is a subject over which there is much 
diversity of opinion. By definition, reflectance factor is re- 
flectance relative to the reflectance of a totally reflecting 
Lambertian reflector, both measured on the same instrument. 
Based on this definition, the best standard for calibrating an 
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Table CII 
Coefficient for the expansion of N g (r)/cosT in powers of the angle of incidence T in radians for three different sphere coatings and three different 
assumptions concerning the self-radiance distribution. 



Coating 



Curve Type 



bo 



b x x 100 



b 2 X 100 



6 3 x 100 



b 4 x 100 



BaS0 4 [14] 



Halon [15] 



BaS0 4 [16] 



upper limit 


1 





7.0 





-1.1 


accepted; 












r<77-/4 
















r>?r/4 


1.03 


-1.240 


-12.4 


16.7 


-5.33 


lower limit 
















upper limit 







5.3 





6.1 


accepted; 












r<7r/4 
















r>7r/4 


1.45 


-219 


405 


-339 


107 


lower limit 
















upper limit 







3.8 





1.6 


accepted; 












r<7r/4 
















r>7r/4 


1.17 


-120 


271 


-248 


79 


lower limit 

















Table GUI 

Values of the integrals in equations (C5), (C6), and (C7) for nine values of 



i 


It 


U 


W 





0.50000 


0.25000 


0.25000 


1 


.39270 


. 12500 


.26770 


2 


.36685 


.07135 


.29550 


3 


.37990 


.04382 


.33608 


4 


.42147 


.02819 


.39328 


5 


.49129 


.01872 


.47257 


6 


.59443 


.01273 


.58170 


7 


.74046 


.00882 


.73164 


8 


.94416 


.00614 


.93802 



Table CIV 
Values of the weighted averages N G (A.) for each of three coatings for three 
wavelengths and three different assumptions concerning the self-radiance 
distribution. 



Determination 
number 



Wave- 
length 



Upper limit 



Accepted 
value 



Lower limit 



BaS0 4 [14] 



Halon [15] 



BaS0 4 [16] 



450 


1.00486 


1.00479 


1.00476 


550 


1.00343 


1.00338 


1.00336 


750 


1.00244 


1.00240 


1.00239 


450 


1.00294 


1.00282 


1.00277 


550 


1.00233 


1.00225 


1.00223 


750 


1.00175 


1.00170 


1.00169 


450 


1.00343 


1.00335 


1.00333 


550 


1.00173 


1.00168 


1.00167 


750 


1.00253 


1.00247 


1.00247 



instrument to measure reflectance factor would be a totally 
reflecting Lambertian reflector, or the closest thing to it 
which could be obtained. The Halon and barium sulfate 
sphere linings which were used in the work described in this 
paper were very close to this ideal. 



However, if reflectance factor measurements are going to 
be of use in practical applications, the instrument must be 
defined as part of the specification of the measurement. It is 
for this reason that we have defined the measurements as 
carefully as possible in terms of integrals describing the 
incident radiation and the instrument response. We believe 
the properties of our reference reflectometer are very close to 
those specified in the definition for directional-hemispherical 
reflectance factor which we have used, and furthermore, we 
believe the geometrical properties of that instrument to be 
stable in time. In order to confirm this, however, we feel that 
it is important to have several working standards representing 
the gamut of the highly reflecting materials which we would 
be measuring with this instrument. It is also important that 
several different types of working standards be used in order 
to be able to detect changes in the optical properties of any 
one of them. As a means of retaining our scale of measure- 
ment in day-to-day measurements and as a means of compar- 
ing one determination with another, we are using the follow- 
ing working Standards: 

Working Standard No. 1 — This standard is a piece of MC- 
20 opal glass purchased in September 1970 from 

Mashpriborintorg 
Smolenskaja pi., 32/34 
121200, Moscow, G-200 
U.S.S.R. 

The piece used is 99 mm by 99 mm and is 20 mm thick. It is 
marked MC-20-2 for identification and is usually referred to 
as the Russian Opal Glass. The side which is measured has a 
very flat, highly polished surface. 

Working Standard No. 2— This standard is a piece of 
Vitrolite glass which was manufactured by 

Libbey Owens Ford 
1701 E. Broadway 
Toledo, Ohio 43605 

The exact date of manufacture of this particular piece of glass 
is unknown, but it has been at least twenty years since any 
glass of this type has been manufactured. The piece used is 
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100 mm by 100 mm and is 11 mm thick. It is marked V6-D1 
for identification. The side which is measured is flat and 
highly polished. 

Working Standard No. 3 — This standard is a porcelain 
enamel on steel plaque which was made around 1946 by 

The Harshaw Chemical Co. 
Division of Kewanee Oil Company 
1933 E. 97th St. 
Cleveland, Ohio 44106 

The plaque used is 108 mm by 108 mm and is marked 80-1 
for identification. The surface which is measured is smooth 
and non-porous, but it has a slight ripple or "orange peel" 
texture. 

These three working standards cover the range of reflec- 
tance factor from nearly 1 to approximately 0.8, which is the 
range generally used for standards to calibrate the gain 
setting of reflectometers. In addition, the Vitrolite standard is 
slightly translucent. Recent tests have shown our instrument 
to be insensitive to this amount of translucence, but if this 
should change the difference will be noticeable in the mea- 



surements obtained using this standard as compared to the 
others. Because of the rippled surface in the porcelain en- 
amel standard, the surface reflection from this standard is 
spread out more than that from the other two standards which 
have flat surfaces. Therefore, differences in the way the 
instrument handles specular reflection may show up in the 
relative measurements of this standard with respect to the 
others. (The most sensitive test of the way the instrument 
treats specular reflectance is made by measuring a specular 
mirror. We intend to use this procedure as soon as we have 
the capability to measure the reflectance of a mirror with the 
necessary accuracy.) 

The cleaning procedure used is an important part of main- 
taining working standards of reflectance. Our present ap- 
proach is to use a procedure which will disturb as little as 
possible the layers of oxides and other surface films charac- 
teristic of the materials. Therefore, we store the standards in 
dessicators and before each measurement we wash them with 
a mild nonfluorescing soap, rinse them thoroughly with hot 
water, rinse them with distilled water, and blot them dry with 
soft tissue paper. Any lint which remains on the surface is 
removed by gently brushing with a soft brush. 
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